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Abstract 


The existence of r]-mesic nuclei in which the rj meson is bound with nucleus 
via the strong interaction was postulated by Haider and Liu over twenty years 
ago, however till now no experiment confirmed it empirically. 

In November 2010, we performed a search for a ^He-r 7 bound state by measur¬ 
ing the excitation function for the dd ^Henvr® and dd ^Hepvr" reactions in 
the vicinity of the rj production threshold. The measurement was performed 
with high statistic and high acceptance with the WAS A detector, installed 
at the cooler synchrotron COSY in the Forschungszentrum Julich. The experi¬ 
ment was carried out using a deuteron COSY beam and deuteron pellet target. 
The beam momentum varied continuously in each of acceleration cycle from 
2.127 GeV/c to 2.422 GeV/c, which corresponds to a range of excess energy 
Q e (-70,30) MeV. 

This dissertation is about the search for ^He-r 7 bound state in dd ^Henyr® 
reaction. The excitation function for the process was determined after identi¬ 
fication of all outgoing particles and the application of the selection conditions 
based on Monte Carlo simulations of r]-mesic helium production and its decay 
via excitation of the N* resonance. The total integrated luminosity was calcu¬ 
lated based on the dd ^Hen and dd pprigpUsp reactions, while the luminos¬ 
ity dependence on the excess energy, used for normalization of the excitation 
function, was determined based on quasi-elastic proton-proton scattering. No 
narrow structure of the r 7 -mesic helium was observed in the excitation function. 
The upper limit of the total cross section for the bound state formation and 
its decay in dd —)• bound —t ^Hen7r° process was determined on the 90% 

confidence level. It varies from 21 to 36 nb for the bound state width ranging 
from 5 MeV to 50 MeV, respectively. However, an indication for a broad state 
was observed. The kinematic region, where we expect the evidence of the sig¬ 
nal from the bound state, cannot be fully described only by the combination 
of the considered background processes. In contrast, the experimental excita¬ 
tion function is very well fitted by the background contributions for the region 
where the signal is not expected. 
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Chapter 1 
Introduction 


For decades, physicists wrestled with basic questions about the surrounding 
universe: What kind of objects it consists of and what kind of interactions 
are responsible for its existence? All matter around us is made of elementary 
particles, which occur in two basic types called quarks and leptons. Unlike lep¬ 
tons, quarks have color charge, which causes the strong interaction. Quantum 
chromodynamics (QCD) is the quantnm field theory describing the strong in¬ 
teractions between quarks and gluons carrying the color charge. According to 
this theory, hadrons consist of three quarks qqq (baryons) or quark-antiquark 
pairs q-q (mesons). The most important baryons are the protons and the 
neutrons, the building blocks of the atomic nuclei. 

One of the most frnitful experimental investigations in the field of nuclear 
physics is the search for new, uncommon objects. Many of them, snch as hy¬ 
pernuclei (^, tetraquarks [^, pentaquarks or dibaryons have been 

already discovered, however still a lot is waiting to be explored. One of those 
theoretically predicted and till now not discovered object is mesic nuclei. This 
new kind of exotic nuclear matter consists of nucleus bound via strong inter¬ 
action with neutral meson snch as r], r]', K, co. One of the most promising 
candidates for such states are the r]-mesic nuclei, postulated by Haider and 
Liu in 1986 [^. The coupled-channel analysis of the ttN —?• ivN , ivN —?• ttttA 
and nN rjN reactions showed that in the close-to-threshold region, the 
r^-nucleon interaction is attractive and strong enongh to form an 77 -nucleus 
bound system |^. However, till now none of experiments confirmed it empiri¬ 
cally. The first theoretical predictions indicated that due to the large number 
of nucleons the 77 meson is more likely to bind to a heavy nucleon, therefore the 
experimental searches concentrated on the heavy nuclei systems. Nevertheless 
those experiments have not brought expected effect . Current researches in- 
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dicate that r} nucleon interaction is considerably stronger than it was expected 
earlier A wide range of possible values of the 'qN scattering length 
calculated for hadronic- and photoproduction of the q meson has not excluded 
the formation of r 7 -nucleus bound states for a light nuclei such as ^He, ^He, 
T and even for deuteron (T^ . 

The existence of mesonic bound state would give unique possibility for bet¬ 
ter understanding the elementary meson-nucleon interaction in nuclear medium 
for low energies. Moreover it would provide information about N*{153b) reso¬ 
nance (Tdl and about q meson properties in nuclear matter [I^. According to 


Bass and Thomas [l6t 17 , the q meson binding inside nuclear matter is very 
sensitive to the singlet component in the quark-gluon wave function of this 
meson, therefore the investigation of the q mesic bound states is important 
also in terms of the understanding of q and q' meson structure. 

It is indicated that a good candidate for experimental search of possible 
binding is ^He-r^ system 
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An observed steep rise in the cross section for 
dd —)■ ^Her^ reaction close to kinematic threshold is a sign of very strong final 
state interaction (FSI), which could be the evidence for the existence of the 
bound system. 

We developed the experimental method which allows for the search for 
^He-r^ bound state in deuteron-deuteron fusion reaction. The proposal for the 
experiment was presented at the meeting of the Program Advisory Commit¬ 
tee in Research Center Jiilich in Germany and accepted for the realization 
in November 2010 (^. The search was performed with high statistic and 
high acceptance at the COSY accelerator by means of the WASA detec¬ 
tion system (To}]^ . The measurement was carried out with deuteron COSY 
beam scattered on internal deuteron pellet target. During each of accelera¬ 
tion cycle the beam momentum was varied continuously from 2.127 GeV/c 
to 2.422 GeV/c crossing the kinematic threshold for the dd —)■ “^Heq reaction 
at 2.336 GeV/c. This range of the beam momenta corresponds to an excess 
energy range from -70 MeV to 30 MeV. The unique ramped beam momentum 
technique allows to reduce the systematic uncertainties. The data were effec¬ 
tively taken for about one week whereof the measurement with magnetic field 
was carried out for only two days because of the failure of cooling system of 
Superconducting Solenoid. 

The search for r]-mesic helium was conducted via the measurement of the 
excitation function for the dd —>■ ^Hen7r° and dd —>■ ^Hep7r“ reactions in the 
vicinity of the q production threshold. The present work is devoted to the 
investigation of the dd —>■ ^Hen7r° reaction. The excitation function for the 
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reaction was determined after the detailed analysis of the experimental data. 
The existence of the bound system should manifest itself as a resonance-like 
structure in the excitation curve for dd bound —t ^Hen7r° reaction 

below the dd —?• ‘^E.er] reaction threshold. In order to interpret the achieved 
experimental excitation functions the advanced Monte Carlo simulations of 
signal dd —)• bound —t ^Hen7r° reaction were carried out. The simulations 

were prepared based on the kinematic model of bound state production and 
decay. According to this model bound nucleus is created in deuteron - 

deuteron collision, rj meson is absorbed on one of the nucleons inside helium 
and may propagate in the nucleus via consecutive excitation of nucleons to the 
A^*(1535) state until the resonance decays into the pion-neutron pair. Before 
the decay, it is assumed that N* resonance moves with Fermi momentum dis¬ 
tribution of nucleons inside ^He. The ^He nucleus, formed from three other 
nucleons, plays then a role of a spectator. The simulations were carried out 
under assumption that the bound state has a Breit-Wigner resonance structure 
with fixed binding energy Bg and a width B and that the beam momentum is 
ramped around threshold for r] production. 

This thesis is divided into ten chapters. The second Chapter presents the¬ 
oretical aspects of search for r^-mesic nuclei. In Chapter 3 the experimental 
background of the search for the r^-mesic nuclei is presented. The fourth Chap¬ 
ter includes general informations about the performed experiment: detection 
facility, the analysis tools, detector calibration and data preselection. The 
Chapter 5 is devoted to the simulations of the dd bound ^Hen7r° 

reaction. Description of the data analysis is presented in Chapter 6 while the 
determination of detection efficiency is presented in the subsequent Chapter. 
Chapter 8 describes the luminosity determination. Chapter 9 presents the final 
results: the excitation function and the upper limit of the total cross section for 
considered process. A summary and the outlook are provided in Chapter 10. 



Chapter 2 

Phenomenology of mesic nuclei 


This chapter is devoted to the overview of theoretical investigations of r^-mesic 
nuclei. The first two sections describe the interaction of rj meson with nucleon 
and the bound states in the scattering theory. In the third section we present 
several predictions for r^-mesic bound states while the fourth section includes 
the physical motivation of the research presented in this thesis. Theoretical 
background including description of the bound and virtual states in scatter¬ 
ing theory basic definitions and formulas are presented in Ref. j^. Detailed 
information reader can also find in the cited literature. 


2.1 r]-N interaction 


The interaction between 77 meson, which properties are presented in Appendix[X| 
and nucleons has been studying since many years paying special attention to 
possibility of the bound states creation. Since, it is impossible to create the rj 
beams due to its short lifetime, the 77 -nucleon studies are based on the investi¬ 
gation of r]N scattering amplitude for the processes like wN rjN, 'yN rjN 
and also NN —>■ NNrj (pp —>■ ppr) j^, pn —>■ pnr] [^). In those reactions 77 
meson interacts with recoiling nucleon and in the low momentum region the 
interaction is dominated by broad nucleon An resonance A^*(1535), which is 
very close to the 77 production threshold (49 MeV above the rjN threshold) and 
has width 150 MeV. The resonance is strongly coupled to the s-wave n — N 
and the 77 — V channels and causes the steep rise in the pion-nucleon cross 
section. Recent and previous experimental data are reviewed in 29 and 
[SOlIM], respectively 


In order to determine the 77 -nucleon scattering amplitude, coupled channel 
calculations have been performed and their results were fitted to the available 
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data. The first calculation carried out by Bhalerao and Liu including rj — N, 
TT — N and A — tt channels results in the strong and attractive interaction be¬ 
tween 77 and nucleon in the low energy (s-wave) region. It was confirmed by 
later calculations |32| - [3^ and allows to postulate possible existence of r^-mesic 
bound states. 


2.2 Bound states in the scattering theory 


The bound state in a usual sense is an object which mass is smaller than the 
sum of its constituent masses. However, in non relativistic quantum mechanics 
binding is more complex. The existence of the unstable states is attribnted 
to the occurrence of poles in the scattering matrix in the complex momentnm 
or energy plane. At the low momenta the scattering matrix can be written 
as (Sbl: 


^ = ( 2 . 1 ) 

1 — ipa 

where p and a are a complex relative r 7 -nucleus momentum and a scattering 

length, respectively. The complex energy E can be expressed by the complex 

2 

momentnm p as E = where is rednced mass of 77 -nucleus system. 


Then the real and imaginary parts are related as Re{E) 


_ Re^ (p)-\-Irin? (p) 

2mu 


and 


Im{E) = 


Re{p)Im{p) 


The pole lying in the physical sheet of momentum and 


energy plane fnlfilling conditions Im{p) > 0 or Re{E) < 0 corresponds to the 


bonnd state or qnasi-bound state, which is schematically presented in Fig. 2.1 


p—plane E—plane 



Figure 2 . 1 : Complex momentum (left) and energy planes (right) with marked 
poles corresponding to bound, quasi-bound, resonant, virtual and quasi-virtual 
states. Figure is adopted from (^. 
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The bound state is related to the case when interaction is described only 
by a real potential (/m(a) = 0). The pole is then located on the positive 
imaginary p axis. In case of inelastic interaction which is associated with 
absorption (Im(a) ^ 0 ) we have the quasi bound state located in the second 
quadrant of the complex momentum plane. The resonance and virtual/quasi- 
virtual state poles he on the unphysical momentum sheet {Im{p) < 0 ) in the 
third and fourth quadrant, respectively. 


2.3 Theoretical predictions for ? 7 -mesic nuclei 


The first theoretical predictions concerning the r^-mesic nuclei existence were 
declared by Haider and Liu in 1986 based on coupled channel calcula¬ 
tions reported by Bhalerao and Liu the year before. Based on the ob¬ 
tained r]N scattering length = (0.28 -|- 0.19i) fm), they postulated the 
formation of r^-mesic nuclei with masses A > 12. However, later phenomeno¬ 
logical and theoretical studies of r] production in hadronic- and photo- in¬ 
duced reactions brought much wider range of possible scattering length from 
(ir]N = (0.18 -|- 0.16z) fm to = (1.03-1- 0.49i) fm 37 . The larger scattering 


lengths do not exclude the formation of a bound states for the helium |ll[| 12 | 
and even deuteron p^ . 

The standard theoretical investigations of the possible binding were focused 
on the construction of the optical potential for the r^-mesic nucleus based on 
information about r]N scattering lengths obtained by fitting the different mod¬ 
els to experimental data and thus, the solution of wave equation. This method 
was used especially in theoretical searches of heavy r^-mesic nuclei using two ap¬ 
proaches j^. In the first approach 77 -nucleon optical potential is constructed 

■—T{r]N r]N)Ap{r 


38 


using "Tp" approximation (Uopt = V + iW = 
where p is reduced pN mass, T{r)N rjN) is p-nucleon transition matrix and 
p is nuclear density). The calculations based on this approach provide informa¬ 
tion about binding energies and widths of p-mesic nuclei for H > 12 [7,39-41 


Another approach is QCD based quark-meson-coupling (QMC), where optical 
potential is constructed with assumption that p is submerged in the nuclear 


medium and couples to quarks and mixes with p' 42,4^. Using this potential 


and solving the Klein-Gordon equation, authors obtained the single particle 
energies for the p meson for closed shell nuclei as well as ®He, ^^B and ^®Mg. 
Obtained results suggest that one should expect bound states in all of those 
nuclei. 

In case of light nuclei, the existence of p-mesic bound states is manifested 
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as poles in the scattering matrix and the corresponding 77 -nuclens scattering 
lengths arf-nudeus- The formation of the 77 -mesic nucleus can proceed if the 
Re{arj-nucieus) is negative, what corresponds to attractive nature of the inter¬ 
action, and the the following inequality is fulfilled |^: 


nucleus^] ^ nucleus^] • 


( 2 . 2 ) 


One of the first predictions concerning light r^-mesic nuclei was carried out 

The author considered rjNN — nNN coupled 


using few body equations 44 


system and observed pole structure corresponding to a quasibound state with 
mass 2430 MeV and width 10-20 MeV. The idea was later used to study possible 
production of d-rj, ^Re-r] and ^He-r] bound states within finite rank approxi¬ 
mation (FRA) 1^. The obtained complex poles in the scattering amplitude 
correspond to the bound states for i?e(o^Ar) G (0.27, 0.98) fm. 

The new approach including information about 77 production mechanism 
and the final state interaction FSI was presented by Neelima Kelkar et al. 

[ 4 ^ . The authors performed analysis of the 77 production, calculated 77 -nucleus 
amplitudes and locate the d- 77 , ^He -77 and ^He -77 mesic nuclei using the concept 
of Wigner’s time delay. This analysis shows, that the formation of light 77 -mesic 
bound states is possible for only small values of while higher scattering 
lengths correspond to resonances |47| . 

Recent phenomenological studies of the ^He -77 bound state production in 
dd —>■ ^Hep 7 r“ reaction were carried out by Wycech and Krzemieh based 
on approximation of the scattering amplitude for two body process. The au¬ 
thors estimated the cross section for dd bound ^Hepyr" process to 

cr ~ 4.5 nb. The result is more than two times higher than the value estimated 


in 1^,48 based on the simple assumption for probability of the hound 

decay in one of possible channel. 


2.4 Motivation for the research 

The discovery of postulated 77 -mesic nuclei would be interesting on its own 
since till now no experiment provides empirical confirmation of its existence. 
The observation of such object would allow to determine its properties and 
thns investigate many important issues in the 77 meson physics. 


One of them are studies of the 77 meson interaction with nucleons inside the 
nuclear matter which would lead to determination of the r]N scattering length 
which is quite poorly known j^, 38 cause cannot be extracted directly 








14 


Phenomenology of mesic nuclei 


from the existing experimental data. 

Moreover, the existence of 77 -mesic nuclei would also provide information 
about A^*(1535) resonance properties in medium m As it was mentioned 
in Sec. | 2 . 1 [ the resonance is coupled to pion and r] meson in the low energy 
region and the bound state studies could provide unique chance to study the 
chiral symmetry of baryons since Ai*(1535) resonance is a chiral partner of 
nucleon 14 . The investigations of 77 -mesic bound states can also be useful in 


testing different approaches related to structure of A^*(1535) resonance [^ , 
cause it is very hard to distinguish between theoretical models from the exist¬ 
ing data 

Another aspect which could be studied via the 77 -mesic nuclei is the struc¬ 
ture of 77 meson. According to [^,43,^ its binding energy is strongly related 
to the contribution of the flavour singlet component of the quark-gluon wave 
function of the 77 meson. The bound states investigation could bring valuable 
information about the magnitude of the glue content in the 77 wave function. 
Moreover, the 77 mass shift inside the nucleus allows to study the axial Ua{1) 
dynamics |4^. 









Chapter 3 

Search for 77-mesic nuclei in 
previous experiments 


The issue of the r^-mesic bound states has become popular already over 25 
years ago when Haider and Liu postulated their existence j^. Since then 
many experiments in different laboratories were dedicated to search for this 
new kind of nuclear matter. The overview of previous measurements carried 
out in heavy and light nuclei regions was described in . This chapter shows 
the summary of experiments and presents current results. 


3.1 Heavy nuclei region 

The first theoretical prediction of the r^-mesic bound states regarded nuclei 
with atomic masses greater than 12 [^[^. Therefore, in the beginning, the 
measurements were performed for the heavy nuclei region. 

First such experiment devoted to the search for r^-mesic nuclei was carried 
out at BNL (Brookhaven National Lab) by measurement oi 7 :^ + A 
p-\-[A— 1)^ reaction with the lithium, carbon, oxygen and aluminium targets. 
Obtained proton spectra did not reveal any peak structure which could be 
interpreted as an indication of the bound state. However, this fact does not 
necessarily mean that the (tt, N) reaction is not a good way to produce r^-mesic 
nuclei. The new investigations with pion beam are going to be performed at 
J-PARC with a new optimal kinematic conditions. It is proposed 

to study of {'K~ ,n) reaction on ^Li and The main advantages over a 

previous BNL experiment will be: (i) back-to-back proton-pion pair detection 
and (ii) the recoilless conditions fulfilled with the pion beam momentum in the 
range between 0.7 and 1.0 GeV/c together with detecting zero-degree neutrons 
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(for the BNL measurement at scattering angle 15° the momentum transfer 
was larger than 200 MeV/c). Moreover, PILOT experiment is planned with 
deuteron target ( 7 r+ + d p + p + rj) in order to estimate background level for 
the considered reactions. 

Another type experiment based on double charge exchange reaction (DCX) 
was performed at LAMPF (Los Alamos Meson Physics Facility) in Los 
Alamos where the r]-mesic ^®F was searched in tt+T —?• 7 r“+ ^®Ne reaction. 

In this case the bound state was considered to be produced via collision of tt'*' 
beam with the neutron inside oxygen nucleus and decays via absorption of the 
r) meson on the neutron and, consequently, the emission of negatively charged 
pion. Obtained excitation functions also did not reveal clear structure which 
could be associated to the r] mesic nuclei. 

The first experiment which claimed an evidence for existence of an r]-mesic 
bound states was performed at LPI (Lebedev Physical Institute) |^[^. The 
r]-mesic nuclei were searched in photoproduction process: 7 +^^C N + {A — 
if) N -\- 71 ^ -\-n-\- X^ where A denotes or nuclei. The invariant mass 
distribution of the correlated 7 r+n pairs shows a narrow peak structure below 
the position of A”* (1535) resonance (shifted by about 90 MeV/c^). The width 
and binding energy of the obtained resonance structure were determined to be 
about 100 MeV and 40 MeV, respectively. Obtained results are in agreement 
with theoretical prediction according to which the production of 77 -mesic nuclei 
proceeds via A* (1535) resonance excitation and its decay into vr-nucleon pair. 

A similar experiment at LPI was dedicated to search for 77 -mesic nuclei 
through observation of the two-nucleon decay mode arising to the two-nucleon 
absorption of the captured 77 in the nucleus [^. In the experiment proton- 
neutron pairs outgoing from carbon target in 7 -l-^^C —^A-l-p-l-n-l-A reac¬ 
tion were measured in coincidence. The protons velocity obtained for photon 
energy £'.y=850 MeV (above 77 photoproduction) peaked at about 0.7 what can 
be interpreted as the result of production of low-energy 77 -mesons followed by 
their two-nucleon annihilation (77 A A —>■ AA). In contrast standard photopro¬ 
duction (for £'.y=650 MeV) does not give the particles with such high momenta. 
Assuming that the observed events from both of described LPI measurements 
(tt+tt and pn decays) are related with the formation and decay of 77 -mesic nu¬ 
clei, the upper limit of the photoproduction cross section was determined and 
is equal to 10 pb. 

The search for back-to-back ■n~p pairs related to the 77 -mesic bound states 
was also carried out at JINR (Joint Institute for Nuclear Research) with the 
internal deuteron beam 1^. The d + ^ 7 r“-|-p-|-A reaction was measured 
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for the deuteron beam energy 2.1 GeV/nucl. In the experiment the 7r~p back- 
to-back correlation was clearly observed and resonance like peak was fonnd 
below the 77 production threshold. The result could be associated with the 
signature of the two-body N* resonance decay related with formation of an 
r^-mesic nucleus. However, the investigation need more intensive beam and the 
higher acceptance of the spectrometer. 

At GSI (Gentre for Heavy Ion Research in Darmstadt) the search for 
r^-mesic nuclei was carried out in recoil-free transfer reaction using similar 
method as in case of measurements of deeply bound pionic states j^. In 
the experiment (d,^He) reaction was measured on ^Li and ^^G targets at GSI 
Fragment Separator System (FRS). The analysis of this data is in progress. 
So far no final result is published. 

A very strong claim for the discovery of the resonance like structure corre¬ 
sponding to the 77 -mesic magnesium was made by the GOSY-GEM group after 
the analysis of p(^^Al,^He) 7 r“p'X reaction |^[^ . Similarly like in case of GSI 
measurement [^, this experiment fulfilled the recoilless kinematics conditions. 
The obtained missing mass spectrum of the ^He shows enhancement for bind¬ 
ing energy of about -13 MeV with the width of about 10 MeV. According to the 
authors, the peak could be interpreted as a signal from ^^Mg-r^ bound state. 
However, it is important to confirm the result with higher statistics. 


3.2 Light nuclei region 


A wide range of possible values of the r^N scattering length a^jv extracted from 
hadronic- and photoproduction of the 77 meson (overview in Ref. |^|^) has 
not excluded the formation of 77 -nucleus bound states for a light nuclei as ^He, 
^He, T (H, 12 and even for deuteron [^. In case of light nuclei 77 absorption 
is smaller and the bound states are expected to be narrower in comparison to 
the case of heavy nuclei. Therefore, the light bound states seems to be good 
candidates for the study of possible binding. 

The experimental studies of the final state interactions (FSI) in ^He 77 and 
^He 77 systems result in observations which may suggest the existence of the 
77 -mesic helium bound states. The measurements performed by SPES-4 [M] , 
SPES -2 g, COSY -11 g and COSY-ANKE g as well as in SPES-4 g, 
SPES-3 jg, GEM [g and COSY-ANKE |g collaborations revealed a strong 
enhancement in the cross section of the dp —>• ^He 77 and dd ^He 77 reactions, 
respectively. This results can be interpreted as a possible indications of the 
77 -mesic helium. Eig. |3.1| shows the cross sections measured for both of consid- 
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ered processes. The fits to the experimental data marked in left and right pan¬ 
els of Fig. M with solid lines gave the value of the 77 -hehum scattering length 
a^3He = [±(2.9 ±0.6) + (3.2±0.4)z] fm and o^4^e = [±(3.1 ± 0.5) + (0.0 ± 
0.5)z] fm (^, respectively. However, these values do not allow to check the 
basic condition for the bound state existence cause it is not possible to verify 
if the real part of the scattering length is larger than the imaginary part. 


D 

§ 400 


0 



excess energy Q [MeV] 


(MeV/c) 


Figure 3.1: (left) Total cross-section for the dp —)■ ^Her^ reaction measured with 
the ANKE (open circles) and the COSY -11 facilities (closed circles) [65] 
and (triangles) j^. Scattering length fit to the ANKE and COSY -11 data 
is represented with dashed and solid lines, respectively, (right) Total cross- 
section for the dd ^Hery reaction as a function of CM momentum obtained 
from the measurements of Frascaria et ah (black diamonds), Willis et 
al. (blue squares), Wronska et ah (magenta triangles) and Budzanowski 
et al. (red circle). The solid line represents a fit in the scattering length 


approximation. The figure is adopted from 62 . 


The COSY -11 and ANKE groups performed additionally measurement of 
the differential cross section for dp —> ^Hery process. The cross section near the 
threshold has not isotropic form because not only s wave but also p wave con¬ 
tributes to the process. It depends linearly on cosd^ and therefore asymmetry 
can be defined as: 


Asymmetry parameter a as a function of ry momentum is presented in Eig. | 3 . 2 [ 
The experimental data were fitted with assumption of very strong variation of 
the s-wave amplitude and not to fast changes of p-wave amplitude . The fit 
is in agreement with measured data and implies the small and constant value 
of the tensor analysing power ^20 for deuteron. The tensor analysing power 
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was recently measured with ANKE group 73 in dp —)■ ^Her^ process. The 


measurement was carried out for the excess energy range Q e (0,11) MeV with 
the polarized deuteron beam. The angle averaged tensor analysing power was 
determined in this region and it varies around -0.2. However, the variation is 
smaller than the error bars what suggest the constant behaviour of ^ 20 - The 
obtained result supports strongly the ESI interpretation in the near-threshold 
region. 



Figure 3.2: The asymmetry parameter a as a function of the p meson mo¬ 
mentum. The data are from (open triangles up), (open triangles 

down), (full squares) and (full dots). The solid curve represents the 
fit according to [^. Figure is adopted from (^ . 


The first direct experimental indication of a light r^-nucleus bound states 
was reported by the TAPS collaboration for the p photoproduction pro¬ 
cess 7 ^He — 7 - TT^pX. The reaction was measured with the TAPS calorime¬ 
ter at the electron accelerator facility Mainz Microtron (MAMI). The mea¬ 
surements of the excitation functions of the 7r°-proton production for two 
ranges of the relative angle between those particles were carried out (upper 
panel of Fig. 3.3). It appeared that a difference between excitation curves for 
opening angles of 170° — 180° and 150° — 170° in the center-of-mass frame 
revealed an enhancement just below the threshold of the y^He —>■ ^Hep re¬ 
action. It was interpreted as a possible signature of a ^He-p bound state 
where p meson captured by one of nucleons inside helium forms an inter¬ 
mediate A^*(1535) resonance which decays into pion-nucleon pair. A binding 
energy and width for the anticipated 77 -mesic bound state were deduced from 
the fit of the Breit-Wigner distribution function to the experimental points 
and are equal to (—4.4 ± 4.2) MeV and (25.6 ± 6.1) MeV, respectively. Those 
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values are consistent with expectations for 77 -mesic nuclei. However, the later 
measurement performed by the TAPS collaboration using upgraded detection 
setup with much higher statistics allows to ascertain that the structure 
observed in the 7 r°-p excitation function is an artefact of the complicated be¬ 
haviour of the background. Obtained results are presented in lower panel 
of Fig. in The excitation functions were measured for the higher photon en¬ 
ergies what allowed to observe the structures corresponding to second and third 
resonance regions of the nucleon. The subtraction of the excitation functions 
for opening angles 165° — 180° and 150° — 165° reveal narrow peak located at 
the 77 production threshold which appears due to the shifting of the low energy 
tails of the second resonance region. 




• 170‘^-ISO'’ 7 

(150*^-170® : 
subiracled) 





Figure 3.3: (upper panel) Excitation functions of the 7 r°-proton production for 
relative angles of 170° — 180° (red triangles) and 150° — 170° (black circles) in 
the 7 °He center-of-mass sytem are shown in the left and center panels. In the 
right panel the difference between both distributions with superimposed line 
denoting the results of the fit of the Breit-Wigner distribution plus background 
are presented. The figure is adopted from [^. (lower panel) Difference be¬ 
tween excitation functions for the opening angle 165° — 180° and 150° — 165° 
obtained by (^. Vertical line corresponds to the 77 production threshold. 
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The cross sections obtained in both analyses 74,7^, presented in Fig. 3.4 


rises very steeply from the prodnction threshold and then stays almost con¬ 
stant. The result is similar to those observed for ^Her^ hadrono-prodnction at 
COSY j^l^. It suggests that the rise of the cross section above threshold is 
independent of the initial channel and is therefore a strong argument for the 
existence of the pole in the scattering matrix which could be associated with 
^He-r^ bound state. 



Fignre 3.4: Total cross section for the y^He —>■ ^Hery reaction. The green 
triangles are from while red dots from . The two vertical lines indicate 
the coherent and the break up thresholds. The inserts show the ratio of data 
and PWIA prediction. Figure is adopted from [76] . 

The search for the ry-mesic ^He was also performed by COSY-11 (77| - |M] 
and COSY-TOF groups via measurement of excitation fnnction of the 
dp pppTT~ and dp —)■ ^He7r° reactions around the rj prodnction thresh¬ 
old. For the first experiment the upper limit of total cross section for dp 

bound ppp7v~ process was estimated to the value of 270 nb and for 
dp —7- bound ^He7r° to the value 70 nb. The analysis of COSY-TOF 

measnrement is still in progress. 

In Jnne 2008 WASA-at-COSY collaboration performed the exclnsive mea¬ 
surement dedicated to search for the ^He-ry bound state in deuteron-deuteron 
fusion reaction. The ry-mesic nuclei was searched via studying of excitation 
fnnction for the dd ^Hep7r“ reaction in the vicinity of ^He-ry threshold. The 
measnrement was carried ont for the beam momentum slowly ramped around 
the ry production threshold corresponding to the range of excess energy Q from 
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about -51 MeV to 22 MeV. Excitation function obtained for the dd —>■ ^Hepyr 
reaction does not show the resonance like structure which could be interpreted 


as a signature of 77 -mesic ^He bound state 25,821. Therefore, an upper limit 


for the cross-section for the bound state formation and decay in the process 
dd — 7 - hound ^Hepyr" was determined at the 90% confidence level. For 

this purpose the excitation function was fitted with Breit-Wigner function with 
fixed binding energy and width combined with second order polynomial. Ob¬ 
tained upper limit presented in Fig. | 3 . 5 | for binding energy 20 MeV varies from 
20 nb to 27 nb as the width of the bound state varies from 5 MeV to 35 MeV. 
The upper limits depend mainly on the width of the bound state and only 
slightly on the binding energy. 

The new data set collected in 2010 with much higher statistics allowed to 
achieve a sensitivity of the cross section of the order of few nb for the bound 
state production in ^Hepvr" reaction. The data analysis for this channel is in 
progress. 


40 


n 

c 


O 

o 

o> 

(0 

E 


0) 

a. 

Q. 

3 


Excluded 



30 35 

r [MeV] 

Figure 3.5: Upper limit at 90% confidence level of the cross section for for¬ 
mation of the ^He -77 bound state and its decay via the dd —> {"^He-r])hound 
^Hepvr” reaction as a function of the width of the bound state. The bind¬ 
ing energy was set to 20 MeV. The green area at the bottom represents the 
systematic uncertainties. Figure is adopted from (^. 










Chapter 4 
Experiment 


This chapter includes general information about the experiment dedicated for 
the search of 77 -mesic ^He which was carried out in 2010. In the first section 
the brief description of the WASA-at-COSY detection system is presented. 
The second section contains information about the tools used in data analysis. 
The last three sections are devoted to accelerator beam settings, calibration of 
appropriate parts of the detector and the data preselection, respectively. 


4.1 Detector Setup 


The experiment described in this thesis was carried out in the Forschungszen- 
trum Jiilich, Germany with the WASA (Wide Angle Shower Apparatus) detec¬ 
tor installed at COSY accelerator. In this section the characteristics of the facil¬ 


ity is briefly described. More detailed description can be found in Ref. |83,84 


4.1.1 cooler SYnchrotron COSY 


The COSY accelerator complex (85| presented in Fig. 4.1 consists of a 184 m 


circumference cooler synchrotron ring, the JULIC cyclotron and the internal 
and external experimental targets. In the COSY ring, protons and deuterons 
(also polarized), pre-accelerated before by JULIC cyclotron, might be accel¬ 
erated in the momentum range between 0.3 GeV/c and 3.7 GeV/c. The ring 
can be filled with up to 10 ^^ unpolarized particles leading to luminosities of 
10^1 
1032 


cm ^s ^ in case of pellet target applied at WASA The beam prepa¬ 
ration includes injection, accumulation and acceleration and takes about few 


cm ^s ^ in case of internal cluster target (ANKE, COSY- 11 ) and 


seconds, while its lifetime with the pellet target (see Sec. 4. 1 . 2 . 1 ) is of the order 
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of several minutes. Beams are cooled by means of electron cooling as well as 
stochastic cooling at injection and high energies, respectively. It allows to 
reach the high beam momentum resolution and decrease the luminosity losses 
during the beam interaction with targets in case of internal experiments. The 
greatest advantage of COSY, in point of view of this work, is the ramped beam 
technique, which permits to perform measurement in slow acceleration mode 
for a given momentum interval within each acceleration cycle (see Sec. 4.3). 
This method allows to reduce the systematic uncertainties which appear in 
case of separate set for each momentum value. 
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Figure 4.1: The COSY accelerator facility with highlighted internal and exter¬ 
nal experimental setups. 
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4.1.2 The WAS A Facility 


The WASA facility is an internal detection system installed at COSY 
since 2007. Before, up to 2005, it was operating at the CELSIUS storage ring 
at The Svedberg Laboratory in Uppsala, Sweden [^. The physics program 
of the WASA-at-COSY facility is dedicated mainly to study of r] and uj rare 
decays [^|^, to the study of dibaryon production and the search for 
r]-mesic nuclei [^|^. The WASA detector vertical cross section is schemati¬ 


cally presented in Fig. 4.2 


COSY 

beam 



The beam and 
target 

interaction point 



FVH 


SEC 


Iron Yoke 


FRH 


Central Detector Forward Detector 


Figure 4.2: Scheme of WASA-at-COSY detection system. The reaction takes 
place in the centre of the detector at the crossing of the COSY beam and 
pellet beam. Gamma quanta, electrons and charged pions being products of 
mesons decays are registered in the Central Detector. Scattered projectiles 
and charged recoil particles like ^He, deuterons and protons are registered in 
the Forward Detector. The abbreviations of the detectors names are explained 
in the text. 


The 47r WASA detector consists of two main parts: Forward Detector (FD) 
and Central Detector (CD) optimized for tagging the recoil particles and reg¬ 
istering the meson decay products, respectively. The internal target of the 
pellet-type is installed in the central part of the detection system (its posi¬ 
tion is marked with dotted line in Fig. 4.2). All individual components of the 
WASA facility are described briefly in the next subsections. 
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4.1.2.1 Pellet Target System 

The Pellet Target system has been developed for the WASA facility to 
fulfil the conditions required for the studies of the rare processes. The main 
components of the system are shown in Fig. |4.3[ 



pellet 
beam tube 

COSY 
beam tube 

pump 


pellet generator 

differential pumping 

scattering chamber 
beryllium tube 
forward window 


pellet 
beam dump 


pump station 


Figure 4.3: The WASA-at-COSY Pellet Target system. 


The Pellet Target setup provides a stream of pellets (frozen droplets) of hy¬ 
drogen (H 2 ) and deuterium (D 2 ). They are produced in the pellet generator, 
located above the Central Detector, where the droplets from the high purity 
liquid jet (H 2 or D 2 ) are formed with the vibrating nozzle. The nozzle vibra¬ 
tions frequency is typically 70 kHz. The droplets freezes by evaporation while 
passing through the chamber becoming the pellets. Then, after the entering 
the 7 cm vacuum-injection capillary the pellets are accelerated up to velocities 
of 60-80 m/s. Finally a skimmer collimates the pellet beam before it enters the 
2 m long pellet tube of 7 mm diameter which is used to guide the pellet beam 
to the interaction region. The average rates of pellets passing the interaction 
point is about few thousand per second. The main pellets characteristics are 
summarized in Table. 14.11 
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pellet size 

~ 35 /rm 

pellet frequency 

5-12 kHz 

pellet velocity 

60-80 m/s 

pellet stream divergence 

~0.04° 

pellet stream diameter at beam 

2-4 mm 

areal target thickness 

> 10^®atoms-cm“^ 


Table 4.1: Pellet Target properties. 


4.1.2.2 Forward Detector (FD) 

The detection and identification of forward scattered projectiles and target- 
recoil particles such as protons, deuterons and helium nuclei and also of neu¬ 
trons and charged pions are carried out with the Forward Detector which covers 
the range of polar angles from 3° to 18°. It consists of fourteen planes of plas¬ 
tic scintillators forming Forward Window Counter (FWC), Forward Trigger 
Hodoscope (FTH), Forward Range Hodoscope (FRH), Forward Range Inter¬ 
leaving Hodoscope (FRI) and Forward Veto Hodoscope (FVH), respectively 
and proportional counter drift tubes called Forward Proportional Chamber 
(FPC). Particles are identified based on measurement of energy loss in the de¬ 
tection layers of FWC, FTH and FRH while their trajectories are reconstructed 
from the signals registered successively in FWC, FPC, FTH and FVH detec¬ 
tors. The registered energy loss permits to determine total particle momentum 
which direction is reconstructed from the measurement of particles tracks by 
means of straw detectors constituting FPC. Components of the Forward De¬ 
tector are presented in Fig. EH and described in text below. 

Forward Window Counter 

The Forward Window Counter (FWC) is the first detector of the Forward Part 
along the beam direction. It consists of two 3 mm layers, each of 24 plastic 
scintillators connected to the photomultipliers (PM) via light guides. The 
FWC is mounted on the paraboloidal stainless steel vacuum window. The first 
layer is shifted by half an element with respect to second one which is mounted 
perpendicularly to the beam direction (see Fig. |4.4(^. The Window Counter 
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is used as a first level of the trigger logic which allows to reduce the background 
coming from the particles scattered downstream beam pipe. It is also one of 
the detector which can be employed in the ^He identification via the AE-AE 
method. 


Forward Proportional Chamber 

The Forward Proportional Chamber (FPC) located directly after FWC is a 
tracking device providing precise angular information about the particles out¬ 
going from the target region (scattering angle resolution about 0.2°). It is 
also used for the accurate reconstruction of the track coordinates of charged 
particles crossing through. The Chamber is composed of 4 modules, each with 
488 proportional drift tubes (straws) of 8 mm diameter made of thin mylar 
foil and filled with argon-ethan gas mixture. The modules are rotated by 45° 
with respect to each other and their orientation is -45°, +45°, 0° and 90° with 
respect to the x direction (see Fig. |4.4(b)|). 


Forward Trigger Hodoscope 

The Forward Trigger Hodoscope (FTH) is third in the order sub-detector con¬ 
sisting of three layers of plastic sintillators. There are 48 radial elements in 
the first layer, closest to the FPC, and 24 elements in the form of archime- 
dian spirals oriented clockwise and counter-clockwise in the last two planes 
(see Fig. |4.4(c) ). The FTH provides for the trigger system angular informa¬ 
tion about the track based on the overlap of hit elements in three consecutive 
layers. Moreover, it gives information about the track multiplicity and is used 
for identification of the charged particles in the FD via energy loss. 


Forward Range Hodoscope 

The five planes of Forward Range Hodoscope are positioned behind the FTH 
(see Fig. |4.4(d)| . Each of them consists of 24 plastic scintillator modules with 
thickness of 11 cm and 15 cm for first three and the last two layers, respectively. 
The energy resolution for particles stopped in the detector is about 3%. The 
FRH together with FWC and FTH allows to determine the energy of charged 
particles stopped in detector or passing through. The initial kinetic energy 
reconstruction and identification of charged particle is based on the energy 
deposited in the different detector planes {AE-AE and AE-E methods). 
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(a) The Forward Window Counter. (b) 3D view of the Forward Pro¬ 
portional Chamber. 



(c) The three layers of the Forward Trigger Hodoscope (left). The 
intersections of elements define pixels as indicated in the projection 
of the planes after hit by two particles (right) 



(d) The Forward Range Hodoscope. 

Figure 4.4: Components of the Forward Detector. 
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Forward Range Intermediate Hodoscope 

The Forward Range Intermediate Hodoscope (FRI) can be situated between 
the second and third layer of the FRH. This two-layer scintillator hodoscope, 
with modules rotated by 90° with respect to each other, delivers precise time 
and position information in two dimensions. This sub-detector was not used 
during the experiment dedicated to this thesis. 

Forward Veto Hodoscope 

The Forward Veto Hodoscope (FVH), being the last subdetector in FD, con¬ 
sists of two layers: one of 12 horizontal and second of 12 vertical plastic scin¬ 
tillator bars. Each bar is equipped with the photomultipliers on both sides. 
The distance between layers is 77 cm. The main goal of FRH is detection of 
high-energetic particles going through the FRH. 

Forward Absorber 

The Forward Absorber (FRA) can be located between the last layer of the 
FRH and the FVH. It is iron plane with thickness of usually 5-10 cm. The 
FRA is used for stopping the slower protons (for example from the pp —)• ppr) 
reaction). The fast protons associated with the elastic scattering penetrate the 
Absorber and induce signals in the FVH which are used for veto purposes in 
the trigger. The absorber was also not used during the described experiment. 

4.1.2.3 Central Detector (CD) 

The Central Detector is built around the interaction point and designed mainly 
for measurements of photons and charged particles originating from 7r° and rj 
mesons decays. The charged particles momenta and reaction vertex are de¬ 
termined by means of Mini Drift Chamber (MDC). Charged particles are 
here bending in the magnetic field provided by surrounding Superconducting 
Solenoid (SCS). First their trajectories are reconstructed, and then knowing 
the magnetic field, the momentum vector is reconstructed. The identification 
of charged particles is based on information about energy deposited by par¬ 
ticles in Plastic Scintillator Barrel (PSB) and in Scintillator Electromagnetic 
Calorimeter (SEC). The calorimeter is also used for the photon identification. 

Mini Drift Chamber 

The Mini Drift Chamber (MDC) is a sub-detector placed around the 60 mm 
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diameter beryllium beam pipe, close to the interaction region and is covered 


by 1 mm thick Al-Be cylinder (see Fig. 4.5(a)). The chamber consists of 1738 


straw tubes arranged in 17 layers covering scattering angles from 24° to 159°. 
The straw diameter is 4 mm, 6 mm and 8 mm in first five inner layers, in 
the 6 middle layers and in 6 outer layers, respectively. The straws are made 
of 25 pm thin aluminized mylar foil and surronnd the 20 pm diameter gold 
plated tungsten anode wire. The first nine inner layers are parallel with respect 
to the beam axis while the next layers are situated with small skew angles 
(6° to 9°). The straws are filled with gas mixture containing argon - ethane 
in ratio 80%-20%. The MDC is immersed in the magnetic field provided by 
the Superconducting Solenoid which causes the bending of charged particles 
trajectories. The main purpose of MDC is determination of particle momenta 
and reaction vertex position. Detailed information about the MDC can be 
found in Ref. |^. 


Plastic Scintillator Barrel 

The Plastic Scintillator Barrel (PSB) is monnted inside the Solenoid Coil and 
surrounds the Drift Chamber (see Fig. 4.5(b)| ). It consists of three parts: 
cylindrical central part (48 scintillator bars) and two endcaps (48 "cake-piece" 
shaped elements each) covering almost dvr solid angle. The main aim of PSB 
is distinction between neutral and charged tracks as well as identification 
of charged particles via AE-E method using total energy information in 
Calorimeter and AE-p method based on momentum information from MDC. 


Superconducting Solenoid 


The Superconducting Solenoid (SCS) installed inside the calorimeter provides 
the magnetic field for the momentum reconstruction of the tracks measured by 
the MDC. The SCS is cooled with the liquid helium and produces the magnetic 


field up to 1.3 T. A detailed description of the solenoid is presented in 94 


Scintillation Electromagnetic Calorimeter 

The Scintillation Electromagnetic Calorimeter (SEC) is situated between the 
SCS and the iron yoke which covers the Central Detector. It is composed of 
1012 sodium-doped Csl scintillating crystals and covers the scattering angles 
from 20° to 169°. The crystals have shape of a truncated pyramid and are orga¬ 
nized in 24 layers. One can distinguish the three main parts of the calorimeter: 
the central with the longest crystals (30 cm), forward build of crystals having 
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25 cm length and the backward consisting of the shortest crystals (20 cm). The 
cross-section of SEC and its angular coverage are presented in Fig. 4.5(c) and 
Fig. |4.5(d) , respectively. The calorimeter is used for the energy measurement 
of charged and neutral particles with resolution of 3% for stopped charged 
particles and 8% for 0.1 GeV photons. Together with MDC and PSB, SEC is 
used for the charged particles identification based on information about their 
deposited energy. Detailed description of the SEC is presented in 



(a) MDC inside the Al-Be cylinder. 


(b) The Plastic Scintillator Barrel. End- 
caps are marked with yellow and red 
colours. 



(c) Cross section of the Scintillating (d) Angular coverage of the SEC. The 
Electromagnetic Calorimeter. numbers above the picture indicate the 

numbers of crystals while their size is 
marked on the vertical axis. 


Figure 4.5: Components of the Central Detector. 
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4.1.3 Data Acquisition System (DAQ) 


The main goal of Data Acquisition system is proper processing of the signals 
from the detector components in order to make them accessible for the data 
analysis. The DAQ system is based on the third generation of the DAQ sys¬ 
tems at COSY and is optimized to cope with the high luminosities (^. The 


structure of WASA DAQ is schematically presented in Fig. 4.6 


event builder 


gigabit fiber 



Figure 4.6: A scheme of the Data Acquisition system (DAQ) for WASA-at- 
COSY. The figure is taken from j^. 

The readout electronic based on Field Programmable Gate Array FPGAs 
used for digitization and buffering of data allows to reach event rates of 10 kHz 
at a live time of 80% of the system j^. Signals from straws and photomulti¬ 
pliers connected with detectors, are distributed and adapted by front-end elec¬ 
tronics (preamplifiers, splitters, discriminators). Subsequently, the analogue 
signals from front-end cards are digitized by means of QDG (Charge-to-Digital 
Gonverter) and TDG (Time-to-Digital Gonverter) read out modules located in 
14 crates. Next, the digitized signals are marked with a time stamps and put 
in FIFO ("First In First Out") queue. Synchronization System (SC), called by 
trigger, computes the event number and send it together with its time stamp 
to all QDC and TDC boards. Signals with a matching time stamp are marked 
with the same event number and passed to the computer readout and to the 
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event builder. Finally, the event are written on the discs. The technical details 


of the DAQ architecture are presented in Refs. (96,97 


4.2 Analysis tools 

For the purpose of investigations made in this thesis, the events generators for 
each of considered reactions were prepared based on proper kinematic models. 
The simulations of the WASA detector response have been carried out with 
the WASA Monte Carlo (WMC) software based on GEANT package (^. The 
data analysis was performed within the RootSorter framework 01 based on 
the data analysis package ROOT developed at CERN (the European 

Organization for Nuclear Research). The ROOT was used for calculations, 
fitting and preparing the histograms shown in this thesis. 


4.3 Beam settings 

The experimental proposal dedicated for the search of bound in 

dd —>■ ^Hen 7 r° and dd —)■ ^Hep 7 r“ reactions with WASA-at-COSY facility was 
accepted for realisation by Programme Advisory Committee in Eorschungszen- 
trum Julich in Germany. The two weeks experiment was carried out at the 
turn of November and December 2010. The data were effectively taken about 
one week, while the rest of time was spent for the accelerator cycle preparation, 
beam and experimental triggers adjustments and pellet target regenerations. 

During the experimental run the momentum of the deuteron beam was var¬ 
ied continuously within each acceleration cycle from 2.127 GeV/c to 2.422 GeV/c, 
crossing the kinematic threshold for the 77 production in the ^Hery reaction at 
2.336 GeV/c. This beam momentum range corresponds to the excess energy 
range of interests Q G (-70,30) MeV ((5 = 0 MeV denotes the threshold). Eor 
the purpose of data analysis this range was divided into 20 intervals. The 
settings of the beam cycle are summarized in Table |4Y , 

The total time of each acceleration cycle in the experimental run has a 
length of 70.3 s. At the beginning of the cycle, the beam is accelerated in 5.7 s 
to the momentum of 2.127 GeV/c via fast ramping. Subsequently, the beam 
momentum is increased slowly from the value of 2.127 GeV/c to 2.422 GeV/c 
and this ramping phase takes 57 s. At time tcyde = 5.5 s the vacuum shutters 
of the Pellet Target are opened and the acquisition system starts recording 
data. In the 63.1’th second of the cycle duration Pellet stream is blocked 
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(shutters are closed) while the data taking continues until 66.1 s. Then the 
DAQ is stopped and the detector voltages are switched off before the beam is 
decreased. 


beam particles 

deuterons 

beam momentum range 

2.127-2.422 GeV/c 

beam cycle length 

70.3 s 

start slow ramping 

5.7 s 

slow ramping time 

57 s 

start DAQ 

5.5 s 

stop DAQ 

66.1 s 

data taking within cycle 

60.6 s (86%) 


Table 4.2: Overview of the accelerator cycle parameters for the experimental 
run. 


4.4 Detector Calibration 


The crucial point in the data analysis of the main considered reactions dd 
^Hen7r° (Chapter]^ and dd —>■ ^Hen (Sec. 8.1) is an identification of ^He ion reg¬ 
istered in the Forward Detector and the determination of its four-momentum. 
The kinetic energy of helium is calculated based on the energy losses in the 
consecutive Forward Range Hodoscope layers. Therefore, it is very important 
to use precise energy calibration of the FRH. The calibration of the FRH (layer 
3 and 4) is described in details in the first subsection. The description of the 
Electromagnetic Calorimeter calibration optimized for the reconstruction of 
photons and hence of 7r° mesons is presented in second subsection. 


4.4.1 Forward Range Hodoscope 

During the data analysis the calibration of plastic scintillator detectors, based 
on the conversion at ADC channels into deposited energy |102 was taken into 
account. However, it was noticed that the calibration is incorrect for Forward 
Range Hodoscope layer 3 and 4 in which high energetic helium outgoing from 


dd —>■ ^Hen reaction is stopped. It is shown in left panel of Fig. 4.7 as dis- 
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agreement in the missing mass spectra of dd ^HeX reaction obtained from 
Monte Carlo simulations and from experimental data. 




Figure 4.7: The missing mass rrix spectrum for dd ^HeX reaction before 
(left panel) and after calibration tuning (right panel). Data is marked with 
red line while Monte Carlo simulation is marked with blue line. 

For the purpose of this analysis, the proper correction for the FRH3 and 
FRH4 calibration was applied. The spectra of energy deposited Edep(FRH3) 
and Edep(FRH4) obtained for data were compared with the spectra obtained 
for WASA Monte Carlo simulations for dd —>■ ^Hen reaction. The comparison 
was carried out for 5 intervals of polar angle in Forward Detector Opo in range 
between 3° and 10.5° and 20 intervals of excess energy Q between -70 MeV and 
30 MeV. Edep(FRH3) and Edep(ERH4) spectra for Monte Carlo simulations 
and data, for each interval, were fitted with gaussian functions in order to 
find the maxima positions - xmc and xp, respectively. Subsequently, peak 
position for data was shifted by offset A to fit the peak position obtained from 
simulations: 


xp — A — Xmc- 


(4.1) 


Exemplary distributions of Edep(ERH3) and Edep(FRH4) for one of the 
chosen 6fd and Q interval with applied fit are presented for Monte Carlo sim¬ 
ulation and data in Eig. |4.8[ Missing mass spectra for Monte Carlo simulation 
and data after all cuts described in Sec. |8.1| with applied calibration tuning fit 
very well, what is shown in right panel of Fig. m 
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Figure 4.8: Spectra of Edep(FRH3) (upper panel) and of Edep(FRH4) (lower 
panel) for data and Monte Carlo simulation, respectively. Eitted gaussian 
functions are marked with red line, while peak positions are marked with blue 
line. The spectra are obtained for 9fd e (3,4.5)° and Q G (-40,-35) MeV. 


4.4.2 Electromagnetic Calorimeter 


The electromagnetic calorimeter is used for measurement of photons emission 
angles and energies. The preliminary SEC calibration was performed based on 
the measurement of cosmic muons and radioactive sources before the WASA 


installation at COSY 95,103 . In order to optimize the photons four-momenta 


reconstruction, the energy calibration was carried out based on determination 
of the invariant mass of neutral pions 7r°. For this purpose events with exactly 
two "neutral" clusters in the Central Detector are selected and considered as 
originating from gamma quanta. Their invariant mass is calculated according 
to below formula: 


^7172 


— cos9i^2)- 


(4.2) 
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where are the measured energies of the photons based on the initial 

calibration while 61^2 is the opening angle between the photons momenta. For 
each gamma quanta pair, being 7r° candidate, the invariant mass is assigned 
to the crystals with the largest energy deposit in the cluster. 

In order to apply a global correction for initial calibration, the distribution 
of invariant mass of two gamma quanta for whole data sample was recon¬ 
structed. The peak position was determined via fitting the sum of signal and 
backgronnd function to the spectrum, what is presented in Fig. |4.9[ 
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Figure 4.9: Invariant mass spectrum for two gamma quanta in Central De¬ 
tector. The red line shows fit to the signal and background while magenta 
line shows polynomial fit to the background. Signal peak after backgronnd 
subtraction is marked as a blue line. The vertical dotted line shows subtracted 
peak position while the vertical solid line shows the 7r° mass. 


The signal was fitted with a Novosibirsk function which is described by 
formula (|4.3|) and (|4.4l). 


P(^x) = / A'^+A'^ 


Qy = 1 + A(x — Xq)/ ax 


sinh{A\/lnA) 

A\/lnA 


(4.3) 

(4.4) 


where: 
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xo-peak position, 
a- width of the peak, 
yl-taih 


The background was fitted as a seven degree polynomial (magenta line). 


The total fit (signal + background) is marked in Fig. 4.9 as a red line, while the 
signal after background subtraction is shown as a blue line. The subtracted 
peak position is depicted as a dotted line and equals 0.1388 GeV/c^. 

The deviation from the actual invariant mass of 7 r° is used to set the values of 
the calibration factor fc, being the ratio of energy for 71 and 72 after applied 
correction to uncorrected energy ( , using the formula: 




m. 


(4.5) 


7172 


The calibration correction factor is applied for each crystal of the calorimeter. 


4.5 Data Preselection 

Data preselection was carried out in two levels: hardware trigger level and the 
selection of the raw data with conditions customized for the present analysis. 
The both levels are discussed in the corresponding sections of this chapter. 


4.5.1 Trigger settings 


The main aim of the hardware trigger system is the reduction of the initial 
event rate to the level that make it possible to be stored on disks, while se¬ 
lecting the events corresponding to the reaction of interest. The trigger con¬ 
ditions are related to multiplicities, coincidences, track matching and energy 
losses in the plastic as well as to the cluster multiplicities and energy deposi¬ 
tion in the electromagnetic calorimeter. In present experiment several triggers 
were set. The main trigger fHedwrl selected events with at least one charged 
particle in the Forward Detector, which corresponds to the track matching 
between FWC, FTH and FRH and in addition with a high energy loss in the 
first layer of the Forward Window Counter. The trigger was dedicated for 
the study of all processes with helium ions in the final state, especially the 


dd —>■ (^He-r^)ftottn(i ^ ^Hen 7 r° (Chapter]^ and dd —>■ ^Hen (Sec. 8 . 1 ) reactions 


The selection of charged particle with a high energy losses, allowed to suppress 
significantly the background coming from fast protons and deuterons in FD, 
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for which the deposited energy is small and close to minimum ionizing particle 
energy loss. 

Additional trigger frhhljpscl was used for the luminosity studies with 
quasi-free pp scattering reaction (Sec. 8.2). It required at least one charged 
particle detected in FD, as well as at least one charged particle in the Central 
Detector. The prescaling factor for this trigger was equal to 1/4000. 


4.5.2 Preselection conditions 

The data selected by the hardware trigger still includes a large sample of 
background events. In order to reduce them and also to limit the computation 
time a preselection of the raw data was carried out. It was performed to select 
only events corresponding to reactions with ^He stopped in Forward Detector 
e. g. dd —>■ bound ^Henyr^ and dd ^Hen. The conditions applied in 

the preselection are following: 

• Exactly one track corresponding to charge particle in Forward Detector, 

• Polar angle of the track 9fd G (3,18)° corresponding to the FD accep¬ 
tance, 

• Graphical cut in Edep(FWCl) vs. Edep(FRHfof) spectrum (energy loss 
in the first layer of Forward Window Counter (FWCl) versus total energy 
deposited in Forward Range Hodoscope (FRH)) in order to reduce the 
background associated with protons and charged pions (Fig. |4.10 ), 

• Based on the Monte Carlo simulations, thresholds for the energy de¬ 
posited in the following layers of Forward Detector were set to the values 
presented in Table |4^ 


FD layer 

threshold [MeV] 

FD layer 

threshold [MeV] 

FWCl 

0.18 

FRHl 

4.0 

FWC2 

0.18 

FRH2 

2.5 

FTHl 

1.5 

FRH3 

2.5 

FTH2 

0.32 

FRH4 

3.5 

FTH3 

0.3 

FRH5 

4.0 


Table 4.3: Thresholds for the energy deposited in the layers of Forward Detec¬ 
tor. 
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Edep(FRHtot) [GeV] 



Edep(FRHtot) [GeV] Edep(FRHtot) [GeV] 


Figure 4.10: Cut applied in Edep(FWCl) vs. Edep(FRHtot) spectrum for: ex¬ 
perimental data (upper panel) and WASA Monte Carlo simulations of the main 
reaction dd —?• i^He-r])hound ^Henvr® (left lower panel), and the normalization 
reaction dd —>■ ^Hen (right lower panel). 




















Chapter 5 

Simulation of the 

dd (4He-77)f,ownd ^ ^HenTr'’ 

reaction 


Present chapter is devoted to the Monte Carlo simulations of the dd 
i^He-ri) hound ^Henvr® process performed based on the kinematic model of the 

r^-mesic helium production and decay. It also includes description of the nu¬ 
cleon momentum distribution inside ^He applied in simulations as well as the 
comparison of these distributions determined for different models. 


5.1 Kinematics of the ry-mesic bound state for¬ 
mation and decay 

We consider the production of the ^He -77 bound state in deuteron-deuteron 
fusion process. The mechanism of the reaction is presented schematically in 
Fig. |5.1[ According to the scheme, the deuteron from the beam hits the deuteron 
in the target. The collision leads to the formation of ^He nucleus bound with 
the 77 meson via strong interaction. The mass of a created bound state is a 
sum of T] and ^He masses reduced by binding energy Bg’. 

bound ~ '^*He ~ Bg- ( 5 - 1 ) 
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DEUTERON 

FUSION 



CREATION OF 
n-MESIC NUCLEUS 


y 



3He, NEUTRON 
AND 27 EMISSION 



ABSORPTION OF n MESON BY 
ONE OF NUCLEON INSIDE THE 
HELIUM 


NUCLEON EXCITATION INSIDE 
THE NUCLEUS- 
N* RESONANCE FORMATION 


Figure 5.1: Scheme of the ^He-r 7 bound state production and decay in 
dd —>■ ^HenTT^ reaction. 


Then, the r] meson can be absorbed by one of the nucleons inside helium and 
may propagate in the nucleus via consecutive excitation of nucleons to the 
A^*(1525) state |104 until the resonance decays into the neutron- 7 r° pair, and 
subsequently 7 r° meson decays into two 7 quanta. It is assumed, that, just 
before the decay, N* resonance momentum distribution can be well approxi¬ 
mated by the Fermi momentum distribution for nucleons inside ^He. The ^He 
plays the role of a spectator which according to the momentum conservation 
in the ^He system moves with the Fermi momentum in the opposite direction 
to the N* resonance. The spectator is considered as a real particle registered 
in the experiment and in the analysis it is assumed that it is on its mass-shell 
during the reaction (iPs^el^ = ^a^e) 105,106 . A very accurate description 


of the r 7 -mesic helium production and decay kinematics including appropriate 
calculations is presented in Ref. m- 


5.2 Simulation scheme 

The simulation of the dd —>■ bound ^Hen 7 r° reaction based on kinemat¬ 

ics presented in previous section, can be schematically described in following 
points: 
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Simulation of the dd —> (^He-ry)bound 


—)■ ^HeriTT® reaction 


1. The deuteron beam momentum Pheam is generated with uniform proba¬ 
bility density distribution in the range of Pheam G (2.127,2.422) GeV/c, 
which corresponds the experimental beam ramping, and then the square 
of invariant mass of the colliding deuterons s^d is calculated from the 
beam and target four-momenta (P^, P^) using the formula: 


sm = |pj + = 2mj (md + \Jmd + IS 


(5.2) 


where: denotes the deuteron mass. 


2. It is assumed that the considered bound state has a resonance-like struc¬ 
ture with fixed binding energy Bg and width B. Therefore, the invariant 
mass of the whole system ^/sBl is distributed randomly according to the 


Breit-Wigner distribution which is given by formula (5.3) and shown in 

Fig-lOl 




rV4 

{\^^dd ~ 'f^{‘^He-r])bound) + -F ^/4 


(5.3) 



Figure 5.2: Breit-Wigner distribution of the invariant mass of the bound 
state system. 
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3. The N* resonance momentum is distributed isotropically in spherical 
coordinates of r^-mesic nucleus {pp, 9*, (J)*) with Fermi momentum dis¬ 
tribution of nucleons inside ^He which is presented for three different 
models in Fig. |5.4| and described in details in next section. 

4. The ^He four-momentum vector is calculated in the CM frame based on 
the momentum conservation and spectator model assumption. 


5. The resonance mass uin* is calculated based on invariant mass y'sdd and 
Fermi momentum pp values according to equation (5.4): 


TfijSl* 


= ( Sdd + -F \p 




(5.4) 


6. The neutron and pion momentum vectors are simulated isotropically in 
the N* frame in spherical coordinates. The absolute value of p**o is fixed 


by the equation (5.5): 


1^** I 

bn,7r0 I 




2 m 


■N* 


(5.5) 


where X{x, y,z) = {x — y — zY — Ayz |108|. 


7. The 7 quanta are simulated isotropically in the frame in spherical co¬ 
ordinates. The absolute value of qs momentum is equal to p*** = mT^of2. 

8. The four-momentum vectors of all ejectiles are transformed into the lab¬ 
oratory frame using the Lorentz transformation. 

9. Simulation of the detection system response is carried out for generated 
events using a GEANT (WASA Monte Carlo) simulation package. 


Fig. 5^ shows spectra obtained for the simulation of the ^He-r] production 
and decay, carried out according to above description, compared with spectra 
related to direct production: dd ^Hen7r°, which is considered as a one of the 
main background contribution in the present studies. 
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^HeriTT® reaction 
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Figure 5.3: The distribution of the 7 r°-neutron opening angle (left panel) and 
the distribution of ^He momentum (right panel), both in CM system obtained 
in simulation of the dd —?• hound ^Henvr® reaction (red line) and of the 

direct dd —>■ ^Hen 7 r° production (black line). The simulation was carried ont 
for the beam momentum range Pbeam G (2.127,2.422) GeV/c. 


For this background process the simulation is performed with an assump¬ 
tion of uniform distribution of the ejectiles over the available phase space. In 
case of dd —>■ {'^He-r])bound ^Hen 7 r° reaction, the relative angle between 
neutron and 7 r° is equal to 180° in the N* reference frame and it is smeared 
due to the Fermi motion by about 30° in the CM frame, while the direct pro¬ 
duction dd ^Hen 7 r° distribution covers the full angular range (see left panel 
of Fig. 5.3). The Fermi motion also determines the ^He momentum in the 
CM system ( 6 *°”J(o and p^^ are strongly correlated). 

In the right panel of Fig. ^ one can see that the ^He momentum distri¬ 
bution in the CM system for the direct reaction is much more wider than for 
the reaction via bound state creation. The cnt in the pf^^ spectrum is used as 
a main criteria in the selection of events corresponding to 77 -mesic helium (see 
next chapter). 

The simulation of dd —?• bound ^ ^Hen 7 r° process was carried out to 

check the feasibility of the measurement of r^-mesic helium production with 
WASA-at-COSY detector setup. These simulations were also used to set ap¬ 
propriate triggering conditions and the beam momentum range dnring the 
experiment, but most of all to compare with experimental data and choose 
the most optimal analysis conditions and cuts (Chapter]^. The Monte Carlo 
simulations have been finally used to estimate the efficiency including all cnts 
applied in the analysis that is presented in Chapter 
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5.3 Nucleon momentum distribution inside ^He 


As it was shown in the previous section, in the simulation of dd —>■ (^He- 
v)bound —^HeuTT^ reaction we assume that N* resonance moves with a Fermi 
momentum given by the distribution for nucleons inside ^He. Unfortunately, 
till now no rigorous calculations for N* momentum distribution inside the 
nuclens are availably The momentum distribution of a particular particle 
depends significantly on the energy which is required to separate this particle 
from the bound system. In ^He there are only nucleons and each nucleon has 
the same separation energy and also the same momentum distribution. In 
case of ^He-r] bound system we can assume that it has binding energy about 
2 MeV, while the cost of nucleon separation is about 20 MeV. If the 77 forms 
an N* resonance with one of the nucleons then the separation energy of the 
N* would be around 22 MeV, and of course the distribution of the N* would 
be very similar to that of a nucleon inside ^He. It works if we assume that 
the mass of the N* is equal to the mass of nucleon and rj meson. In fact 
the average mass of the N* in vacuum is much higher which would imply a 
rather different distribution peaked at higher values. However, in reality we 
don’t know the mass of the N* inside a nucleus and it is difficult to find good 
solution [lo^ . According to [IIol[m] N* resonance momentum distribution 
inside the ^He-r^ bound state can be in a good approximation described by the 
momentum distribution of neutron inside ^He and therefore we apply it in our 
simulations. 

The Fermi momentum distributions for nucleons inside atomic nuclei are 
calculated based on different interaction models. For nucleons inside ^He the 
momentum distributions predicted by three independent models are shown in 
Fig. |5.4[ The distribution represented by a thick line is calcnlated from helium 
wave function derived based on Fermi three parameter charge distribution of 
nucleus |25tpT^. The momentum distribution is described by the formula 


(5.6): 


fip) = -exp(^ 
a \ b 


(5.6) 


where a = 0.0001989184519 (GeV/c)^ b = 0.0028615450879 (GeV/c)T Fermi 
momentum is given in units of GeV/c. 


first theoretical calculations are at present carried out by N. Kelkar 46 
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Simulation of the dd —> (^He-ry)bound 


—)■ ^HeriTT® reaction 



Figure 5.4: Fermi momentum distribution of nucleons inside ^He given by 
analytic formula (thick solid) and estimated for the AV18-TM model (dashed) 
and the CDB2000-UIX model (dotted). The distributions were normalized to 


unity in the momentum range from 0 to 0.5 GeV/c 107 


The formula is determined based on 3-parameter charge distribution of nu¬ 
cleus |112| : 


f r - c 
1 -|- exp - 

V ^ 

with parameters: 
w=0.445±0.020 
c=(1.008±0.013) fm 
2=(0.327±0.002) fm. 


p(r) = po 1 + 


wr 


(5.7) 


The dashed and dotted lines depict distributions obtained from AV18 and 
the CDB2000 nucleon-nucleon interaction models in conjunction with Ur- 
bana IX (UIX) and Tucson-Melbourne (TM) three nucleon interaction (TNI) m 
Due to the fact that ^He is symmetrical nucleus, proton and neutron momen¬ 
tum distributions are in good approximation equal. 

In the low momentum region pp up to 0.4 GeV/c it is no difference between 
distributions determined using 3 Nuclear Force’s 3NF’s like AV18-TM and 
CDB2000-UIX (Fig. |5.4[ ). A slight discrepancy is visible for pp >0.4 GeV/c 
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and results from different interaction Hamiltonian forms defined for above- 
cited models. However, the difference between the distributions derived from 
AV18 and CDB2000 models and given by analytic formula is significant and 
the maxima of these distributions are shifted by about 45 MeV/c |107| . The 
discrepancy can result from the fact that the formula (5.6) was derived from 
nucleus charge distribution smeared out by the charge distribution of protons, 
whereas the AV18 and the CDB2000 models allow for the finite size of nucleus 
charge distributions and are related to the momentum of the point like protons 
in the alpha particle [IIo| . 


The simulation of dd bound ^Hen7r° —>■ ^Henyy reaction was 

carried out for each of the above mentioned models according to description in 
Sec. |5.2[ The geometrical acceptance of WAS A detector as a function of excess 
energy Q for the considered models is presented in Fig. |5.5[ The acceptance 
was determined for simultaneous registration of ^He in Forward Detector and 
two y quanta in Central Detector. 



Figure 5.5: Geometrical acceptance of the WASA-at-COSY detector in case 
of dd —)• bound ^Hen7r° ^Henyy reaction for the three different 

models of nucleon Fermi momentum distribution inside ^He: analytic formula 
(black line) AV18-TM (red line) and CDB2000-UIX (blue line). The simulation 
was carried out for the bound state width r=2h MeV and binding energy 
B,=10 MeV. 
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Simulation of the dd —> (^He-ry)bound 


—)■ ^HeriTT® reaction 


We can see significant difference between result obtained for analytic for¬ 
mula (5.6) and the models AV18-TM and CDB2000-UIX. The acceptance de¬ 
termined from the simulation with the analytic formula is higher of about 15% 
than the acceptances for the other two models. It follows that usage the dif¬ 
ferent models describing the nucleon momentum inside ^He gives one of the 
highest contribution to the systematic errors in the data analysis. However, it 
is important to stress that the shape of the acceptance is model independent 
and therefore the condition about possible existence of the ^He-r] mesic nuclei 
does not depend on the model of the Fermi momentum distribution. In our 
studies we will consider more realistic nucleon momentum distributions calcu¬ 
lated from AV18 and the CDB2000 potentials. 




Chapter 6 

Analysis of the ^HenTr*^ events 


The analysed data set consists of 810 runs (run numbers: 22758-22889, 22994- 
23025, 23254-27399) and corresponds to an effective measurement time of 
about 155 hours. Most (about 84%) of the measurement was carried out with¬ 
out magnetic field in CD due to the failure of the Solenoid cooling system. The 
main part of the analysis was devoted to selection and reconstruction of events 
corresponding to the dd —)• C^He-r])bound ^Hen7r° process. This analysis is 
described step by step in the following sections. 


6.1 Events Selection 


The initial events selection for the considered reaction is carried out on the 
hardware trigger level which requested at least one charged particle in the FD 


and a high energy deposition in the FWC detector (Sec. 4.5.1). Subsequently, 
the preselection dedicated for ^He is performed in order to speed up the anal¬ 
ysis, as it was described in Sec. |4.5.2[ In the next step of the analysis, all 
ejectiles are identified and the events, which may correspond to the produc¬ 
tion of bound states, are selected with appropriate cuts based on the Monte 
Carlo simulations. ^He is registered in the Forward Detector, while gamma 
quanta from the 7r° decay in the Central Detector. Angular distributions for 


the outgoing ^He and q’s are shown in Fig. 6T. An angular ranges covered by 
respective parts of WASA-at-COSY detection setup are marked with shaded 
areas. The scheme of WASA detection setup with marked dd —>■ {'^He-r])bound 


^Hen7r° process is presented in Fig. 6^, The methods of particles identification 
are described in next subsections. 
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e[deg ] 



0[deg ] 


Figure 6.1: Simulated angular distributions of ^He (a), and gamma quanta 
(b) outgoing from the dd (^He-r7)bs —?• ^Henvr® ^ ^Henqq reaction. Figure 
shows results for Monte Carlo simulation generated using the AV18 potential 
model describing momentum distribution of nucleons inside ^He. 
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Figure 6.2: Scheme of the WASA-at-COSY detection system with tagged dd 
i^He-ri)hound ^Henvr® reaction. Helium is registered in Forward Detector 

whereas gamma quanta are detected in the Central Detector. 


6.1.1 ^He identification in the Forward Detector 

According to performed simulations of dd —)• {'^He-r])bound ^Hen7r° reaction, 
^He is mostly (95%) stopped in first layer of Range Hodoscope in Forward 
Detector and just in 5% in the second layer. It is presented in left panel of 
Fig. O where the spectrum of ^He momentum in the CM is plotted for helium 
stopped in FRHl and FRH2 layers. 

In order to reduce significantly background originating from higher ener¬ 
getic helium (right panel of Fig. 6.3) with just small signal reduction, the veto 
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p-JGeV/c] 



Figure 6.3: Momentum distribution of ^He in the CM system obtained from 
simulations (left panel) and experimental data (right panel). The red, blue and 
black lines denote momentum distribution of all registered ^He and stopped in 
FRHl, FRH2, respectively. 


condition was set on the second and further FRH layers (deposited energy 
Edep(FRH2,3,4,5)<0.015 GeV). The ^He was identified with AE-AE method 
based on energy losses in the FWCl and FRHl. In Fig. |6.4| one can see the 
spectrum of the Edep(FWCl) vs. Edep(FRHl) with marked graphical cut 
applied for ^He ions selection. 



^ 0.02 
® 0.018 
“0.016 
O 0.014 

i 0.012 

a 0.01 
S 0.008 
0.006 
0.004 
0.002 

° 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Edep(FRH1)[GeV] 


900 

800 

700 

600 

500 

400 

300 

200 

100 

0 


Figure 6.4: Experimental spectrum of energy deposited in FWCl and FRHl. 
The selected area for ^He is marked with black line. The empty area below 
comes from the preselection cut (see Sec. 4.5.2). 
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6.1.2 7r° and neutron identification in the Central Detec¬ 

tor 


As it was mentioned in Sec. 6.1, 7r° mesons from dd —?• hound ^Hen7r° 

reaction are registered in the Central Detector. The neutral pions 7r° are 
reconstructed from the invariant mass of two gamma quanta originating from 
its decay. In the analysis, hrst all events with at least two neutral clusters in 
electromagnetic calorimeter are selected. Next, all gamma pair combinations 
are considered, and for each pair the invariant mass rrij^ is calculated. In 
case of more than two clusters, we take into account only this combination of 
clusters for which the difference between mass m^rO and invariant mass of two 
gamma quanta is minimal. The cut applied in invariant mass spectrum, 
based on Monte Carlo simulations, is presented in Fig. |6.5[ Experimental data 
is marked with the red line, while the Monte Carlo simulations with the black 
line. 


— DATA 
—Simulations: 





0.4 

rriy,^ [GeV/c^] 

Figure 6.5: 7r° identification via cut in invariant mass spectrum. Applied cuts 
are marked with green lines. 


Neutron was identified via the missing mass technique. Knowing a four- 
momenta of deuteron beam (E^, ph), deuteron target (Ej = rrid, Pt = 0), 
helium {Esjje, Pme) and 7r° (E^o, p'^o) and employing the principle of momen¬ 
tum and energy conservation we can calculate the missing mass as follows: 


= Ex-Px = {Ed + E] - E3He - E^of - {Pb +Pt- psHe “ pZ^f- (6-1) 

Unfortunately the missing mass spectrum contains a lot of background 
from reactions with more than two gamma quanta in the decay channel (most 
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probably with the dd —)■ ^Hen 7 r° 7 r° reaction). This background has been sig¬ 
nificantly reduced by applying the momentum cut on the calorimeter clusters 
(third or more), which were not selected as 7 coming from 7 r° decay, by the 
reconstruction procedure. The cut is based on Monte Carlo simulations, and 
for the further analysis only these events are accepted for which momentum 


corresponds to additional cluster is less than 0.03 GeV/c as shown in Fig. 6.6 




Figure 6 . 6 : Momentum distribution for additional neutral clusters for Monte 
Carlo simulation (left panel) and data (right panel). The applied cut is marked 
by red line. 


Additionally, the cut on the mx{Ex) spectrum was applied as shown in Fig. 6.7 
to reduce the background coming from dd —>■ ^Hen 7 r° 7 r° reaction. 



Figure 6.7: Missing mass rrix vs. missing energy for Monte Carlo simulation 
(left panel) and experimental data (right panel). The applied cut is marked in 
black. 
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The plot comparing missing mass spectra for data before and after back¬ 
ground reduction cuts is presented in the left panel of Fig. 0 while the 
neutron spectrum after final cuts is shown in the right panel. Most of the 
remaining background at the right side of the spectrum was rejected applying 
cuts marked with vertical lines. 




Figure 6.8: The missing mass spectra for dd —?• ^HeX reaction, (left) Com¬ 
parison of the missing mass nix before and after applied cuts for experimental 
data, (right) The final missing mass spectrum. Data is marked with red line, 
the Monte Carlo simulations of the signal is marked with the black line. The 
region accepted for further analysis is marked with vertical lines. 


6.1.3 Kinematic cuts for (^He-r7)&ound events selection 


As it was mentioned in Sec. 5.1, the ^He in the dd —>• bound ^Hen7r° 

reaction plays a role of the spectator which is moving with the low momentum 
corresponding to the Fermi momentum of the nucleons inside ^He. Therefore, 
in the spectrum we selected two regions: region where we expect a sig¬ 
nificant contribution form the bound state signal for G (0.07,0.2) GeV/c 
and the region poor in signal where the background dd —>■ ^Hen7r° and dd 
^HeA* ^ ^Hen7r° processes are dominating (region G (0.3,0.4) GeV/c). 
These regions referred to as "Signal Rich" and "Signal Poor" are marked in 
the left upper panel of Fig. |6.9| as a region A and B, respectively. 


In order to improve the selection of events corresponding to the (^He-p) 


bound 1 


an additional cuts reducing the background contributions, were applied in the 
neutron and pion kinetic energies in the CM system. The energy spectra with 


marked cuts are presented in the upper right and lower left panels of Fig. |6.9 
Based on the spectrum obtained from simulations of dd —>■ bound 

^Hen7r° process and presented in the left panel of Fig. 5.3, we applied also a 
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cut in the neutron- tt® opening angle in the CM frame corresponding to the 
range between 145° and 180°. Since the opening angle is strongly correlated 
with the ^He momentum the cut removes only a small amount of events below 


145° what is visible in the right lower panel of Fig. 6.9 



pf„";[GeV/c] EKi„„[GeV] 



Figure 6.9: Spectrum (left upper panel), distribution in region A 

(right upper panel), distribution in region A (left lower panel) and o 

distribution in region A (right lower panel). Data are shown in red. Monte 
Carlo simnlations of signal are shown in black, while the applied cnts are 
marked with the green lines. 





















Chapter 7 

Detection efficiency 


The experimental data are collected with non perfect geometrical acceptance as 
well as non-perfect detection and reconstruction efficiency. In order to correct 
the obtained results for those detector effects, one should carefully study the 
behaviour of acceptance and efficiency distributions. 

The overall detection and reconstruction efficiency was determined based 
on the Monte Carlo simulation for the dd —>■ {"^He-r])bound —t ^Hen7r° process 
carried out taking into account detection system response and all selection cuts 
described in Chapter The efficiency was calculated as a ratio of the number 
of events accepted by detection system to the number of generated events. 
The correction of the experimental data is applied by dividing the determined 
distributions of observables of interests by the full efficiency. The efficiency 
for the "Signal Rich" region A (see Sec. 6.1.3), together with the detector 
acceptance are presented as a function of the excess energy in Fig. |7.1[ It is 
worth to emphasize that the efficiency does not depend on the bound state 


width r and the binding energy Bg as it is shown in Fig. 7.2 


The geometrical acceptance of the WASA-at-COSY detector for the dd 
i^He-ri)bound ^Henvr® reaction is equal to about 57% while the full efficiency 

including all cuts applied in the analysis is about 9% and is smooth in the 
whole excess energy range. 
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lU 31 

Q [MeV] 


Figure 7.1: The acceptance and efficiency for the registration and reconstruc¬ 
tion of dd —7- bound ^ ^Henvr® reaction as a function of excess energy Q. 

The geometrical acceptance of the WASA detector is shown with red triangles 
while the full efficiency including detection and reconstruction efficiency for 
the region rich in signal is shown with black circles. 



Figure 7.2: The efficiency for the registration and reconstruction of dd —?• (^He- 
v)bound ^Hen7r° reaction for the region rich in signal as a function of excess 
energy Q. Results obtained for different bound state widths F and the binding 
energies Bg are marked with different colours. 














Chapter 8 


Luminosity Determination 


In this chapter two methods of the luminosity determination are presentecQ 
As it was described in Sec. |4.31 the technique of continuous change of the beam 
momentum in one accelerator cycle was applied in the experiment. During an 
acceleration process the luminosity could vary due to beam losses caused by 
the interaction with the target and with the rest gas in the accelerator beam 
line, as well as due to the changes in the beam-target overlap correlated with 
momentum variation and adiabatic shrinking of the beam size 114 . There¬ 
fore, it is necessary to determine not only the total integrated luminosity but 
also its dependence on the excess energy. The total integrated luminosity is 
determined based on the dd ^Hen and quasi free pp —>■ pp reactions for 
which the cross sections were already experimentally established. Because of 
the acceptance variation for the beam momentum range for which ^He ions are 
stopped between two Forward Detector layers, the excess energy dependence 
of the luminosity is determined based on quasi-free pp —?• pp reaction for which 
the WASA acceptance is a smooth function of the beam momentum. 

The precise luminosity determination as a function of excess energy Q 
is important for the normalization of the obtained excitation function for 
dd — 7 - ^Hen 7 r° reaction and hence for the interpretation of the result in view of 
the hypothesis of the ^He -77 bound state production. 


^The description of luminosity determination has been already published by the author 
in a form of conference proceedings (Acta Phys. Polon. B46 (2015) 1, 133) 
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8.1 Integrated luminosity - dd ^Hen reaction 
analysis 


Cross section determination 

The absolute value of the integrated luminosity was determined using the ex¬ 
perimental data on the dd —?• ^Hen cross-sections measured by SATURNE col¬ 
laboration for four beam momenta in the range between 1.65 and 2.49 GeV/c ini 
The cross section add^^Hen dependence on the square of the momentum trans¬ 
fer t = — Pfoeam)^ niaj be parametrized as follows fTT^[Tl6| : 


) , ( 8 . 1 ) 

1=1 

where t^ax denotes maximal value of t measured for a given beam momentum 
at SATURNE. Parameters aj and bi are described as a function of the total 
energy 


pciri{y/^d) = 


Pi 


+ Ti 


( 8 . 2 ) 


Qi 

where the values of Pi, and fj were determined in by the fit of the above 
formula to the cross sections measured at SATURNE n. The cross section 
parametrization was described in details in mi The parameters obtained 
from the fit to SATURNE data |115| are shown in the Table 8.1, 



Pi 

Qi 

n 

ai 

11.64 

4.05 

-14.49 

hi 

0.78 

3.92 

9.04 

02 

2327.04 

-1.44 

-399.27 

b2 

0.78 

3.92 

9.04 

as 

0.22 

4.08 

1.24 


0.78 

3.92 

9.04 


Table 8.1: Parameters for the t and ^/^ dependence of the total cross section 
of the dd ^Hen reaction. In the applied parametrization, the t and 
values are expressed in (GeV/c)^, GeV and pb/(GeV/c)^, respectively. 


The differential cross section as a function of t — tmax for the three different 
beam momentum values from our experimental range p^eam G (2.127,2.422) GeV/c 
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and the total cross section as a function of the invariant mass are presented in 
the left and right panel of Fig. in respectively. 



t-t„3, [(GeV/c)^] [GeV] 


Figure 8.1: (left) Differential cross section for p^eam = 2.127, 2.334, 2.422 GeV/c 
and (right) total cross section as a function of the - dashed red line covers 
the experimental beam momentum range Pbeam G (2.127,2.422) GeV/c while 
green line shows the threshold for p meson production. 


We may determine angular dependence of the cross section using a following 
relation: 


da da dt 

d{cosO*) dt d{cos6*) 

with the Jacobian term = 2 • • IPa^el calculated based on the 

momentum transfer squared in the GM system: 

t = {Fsffg — FbeamY = ^ • F'bearn ' ^ ' \Pbeam I ' \P^He I ' > (8-4) 

where ^mei Pheami Pme beam and ^He energy, momenta and 

the ^He emission angle in the CM frame, respectively. 

The relation between the ^He scattering angle 6iab and cosd* is presented 
in Fig. |8.2[ The ^He angular range from about 4° to 10° corresponds to the 
cose* e (0.88,0.98). 

The available SATURNE experimental data closest to the range of beam 
momentum used in the experiment for the angular range relevant for our anal¬ 
ysis are shown in the left panel of Fig. |8.3[ Superimposed lines present re¬ 
sults of the above described parametrisation for beam momenta corresponding 
to the experimental points: 1.992 GeV/c and 2.492 GeV/c (red and black, 
respectively) and for two exemplary momenta corresponding to Q = 0 and 
Q = -40 MeV. 
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0iab [deg] 


Figure 8.2: cosO* dependence on the 9iab angle for dd —)■ ^Hen reaction (red 
line). The horizontal lines mark the selected intervals of cosO*. 




Figure 8.3: (left) Differential cross section as a function of cos9* for SAT- 
URNE experimental data (squares/red and dots/black points for fixed beam 
momentum Pbearn=l-992 GeV/c and Pbeam=2.492 GeV/c, respectively) and ob¬ 
tained from parametrization (top solid/red, dashed/green, dash-dotted/blue, 
and bottom solid/black lines for Pbeam equal to 1.992 GeV/c, 2.217 GeV/c, 
2.335 GeV/c and 2.492 GeV/c, respectively), (right) The ratio of experimen¬ 
tal and parametrized cross section —! d(cos 0 *) P6earn=1-992 GeV/c 

(squares/red) and PbeaTO=2.492 GeV/c (dots/black) fitted with second degree 
polynomial functions (dashed/red and solid/black lines, respectively). The 
marked errors result from the statistical experimental nncertainties. 


In the angular region of interest the experimental points lie below the 
curves obtained based on the parametrization defined in Eq. (8.1) and (8.2). 
The discrepancy significantly affects the luminosity determination, therefore 
correction for the parametrization was necessary and was applied for the 
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CO50* G (0.88,0.98). The ratio between experimental and parametrized cross 
section fitted with a second degree polynomial function 

for both experimental beam momentum values: 1.992 GeV/c and 2.492 GeV/c. Ob¬ 
tained result is presented in the right panel of Fig. |8.3[ The cross section correc¬ 
tion is calculated for fixed using the fitted functions and linearly interpo¬ 
lated for the proper beam momentum value from range ^ (2.127,2.422) GeV/c. 


Selection of dd ^Hen events 

The measurement of the dd ^Hen reaction was based on the registration of 
the outgoing helium in the Forward Detector. In the first step of analysis at 
least one charged particle in the FD and a high energy deposition in the FWG 
Detector (Sec. 4.5.1) were required to reduce the background especially from 
protons and pious. Then the preselected data (see Sec. 4.5.2) were taken into 
account. 
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Figure 8.4: (upper panel) Edep(FRHl) vs. Edep(FRH2) spectrum, (left lower 
panel) Edep(FRH2) vs. Edep(FRH3), (right lower panel) Edep(FRH3) vs. 
Edep(FRH4). The spectra correspond to events for ^He stopped in FRH3 or 
FRH4. Solid line shows graphical cut used in the analysis for the selection of 
events with ^He ions. 
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Low-energetic ^He ions were stopped in the 3rd layer of the Forward Range 
Hodoscope, while high-energetic in the 4th layer. The helium identification 
was based on the AE-AE method as presented in Fig. |8.4[ In order to disen¬ 
tangle the ^He from other charged particles in FD, a cut in the Edep(FRHl) 
vs. Edep(FRH2) spectrum was applied (upper panel of Eig. |8.4| ). Next, he¬ 
lium stopped in ERH3 or in ERH4 was selected with cuts in Edep(FRH2) vs. 
Edep(ERH3) and Edep(FRH3) vs. Edep(ERH4), presented in left and right 
lower panels of Fig. |8.4[ respectively. 

The outgoing neutrons were identified using the missing mass technique. 
In order to reduce the background originating from the multi-pion reactions 
like dd ^Hen7r°7r° the number of neutral clusters reconstructed in CD was 
requested to be less than 2. Then, to reduce background contribution arising 
from quasi-free dpijisp) ^He(nsp)7r°, the cut in missing mass nix vs. missing 


energy E^ spectrum was applied as it is presented in Fig. 8.5 



Ex [GeV] 




ExIGeV] 


Figure 8.5: Spectra of missing mass rrix vs. missing energy E^ for simu¬ 
lation of dd —>■ ^Hen reaction (left upper panel), simulation of background 
dpijisp) ^He(nsp)7r° reaction (right upper panel) and data (lower panel). Ap¬ 
plied cut is marked with black line. 
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Additional, for the high beam momentum region background was sub¬ 
tracted via fitting the signal and background function to the missing mass 
spectrum for different intervals of cos9* and beam momentum. An example of 


the result for Q E (0,5) MeV and cos9* E (0.96,0.98) is presented in Fig. 8.6 


The signal is described with a Novosibirsk function, which is given by for¬ 
mula (4.3) and (4.4) (see Sec. 4.4.2). The background is fitted with a Gauss 
function. Fit to the signal and background is marked as a red line, while the 
background alone is marked as a green line. 
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Figure 8.6: The missing mass spectrum for cos9* E (0.96,0.98) and 
Q E (0,5) MeV. The red line shows fit to the signal and background while 
green line shows fit of the Gauss function to the background. Signal peak is 
marked as a blue line. 

For the low beam momentum regions, the obtained missing mass spectra 
were almost background-free, therefore no fitting procedure was necessary. 

Luminosity determination 

In order to calculate the total integrated luminosity, the number of events N, 
the efficiency eij, as well as cross section was determined for 5 intervals 

of cos9* in the range from 0.88 to 0.98 and 5 intervals of excess energy Q in the 
range from -70 MeV to 30 MeV. The integrated luminosity was then calculated 
for each (f, j)-th interval in following way: 

N- ■ 

j^int _ _ ^ _ 


(8.5) 
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where A{cos9*) is the width of the cos9* interval. The overall efficiency in¬ 
cluding reconstruction efficiency and geometrical acceptance of the detector 
was determined based on the Monte Carlo simulations and it varies between 
50% and 70%. 

The luminosity dependence of cos9* for whole excess energy range is pre¬ 
sented in Fig. |8.7[ 



Figure 8.7: Integrated luminosity as a function of cos9*. The statistical un¬ 
certainties are marked as vertical bars. The weighted average of integrated 
luminosity is marked as a dashed red line and is equal to 1102 ± 2 nb“^ where 
only a statistical error is given. 

The total integrated luminosity was calculated as a weighted average of the 
luminosities determined for individual cos9* intervals: 


T tot _ 

^ dd^^Hen 


eLi u 


ALiY 




{ALiY 


AL 


tot _ 

dd^^Hen 



- 1/2 


( 8 . 6 ) 


The average integrated luminosity with its statistical uncertainty is equal to 


L'dLmen = (1102 ± 2) nb-1 (see Fig. 8.7). 


8.2 Dependence on the excess energy - quasi- 
free dd ^ ppnspTisp reaction analysis 


In order to determine the luminosity dependence on the beam momentum 
we used the quasi-elastic proton-proton scattering in the deuteron-deuteron 
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collisions: dd PpngpUgp, which, in contrast to dd ^Hen reaction, is char¬ 
acterized by the smooth acceptance in the whole momentum range. In this 
reaction protons from the deuteron beam are scattered on the protons in the 
deuteron target. We assume that the neutrons are acting only as spectators 
which means that they do not take part in reactions. 

In the case of quasi-free proton-proton scattering the formula for the calcu¬ 
lation of the integrated luminosity can be written in the following form (ml: 

_ ^O^exp /Q yN, 

^dd^ppUspUsp 27 tI ^ ) 

where: 


/ = r — 

J AQ{Blab,<l>lab) dQ 


{0 *, (j )*, PFf ^, Opf ^, (pFb ) PFt ) (^Ft ) 4>Ft ) 


^ f ipFbi ^Fb) 'pFbiPFtt (^Ffi (l)Ft)dpFbdcos6FbdcjyFbdpFfdcosOFbd<f)Fbd(l)*dcosO*. 


The formula is determined based on the fact, that the number of quasi-free 
scattered protons into the solid angle AQ{9iab, (piab) is proportional to the inte¬ 
grated luminosity L, as well as the inner product of the differential cross section 
for scattering into the solid angle around 6* and <p* expressed in proton-proton 
CM system: ^ {6 *, 0*, pp ^, Opb , 0^6; PFt > (^Ft: 4>Ft ) and the probability density of 
the Fermi momentum distributions: • f{pFt,OFt,(t>Ft) inside the 

deuteron beam (6) and deuteron target (t), respectively. The detailed de¬ 
scription of the luminosity calculation for quasi-free reaction one can find in 
Ref. firn . 

Due to the complex detection geometry a solid angle corresponding to par¬ 
ticular part of the detector cannot be in general expressed in a closed analytical 
form. Therefore, the integral in above equation was computed with the Monte 
Carlo simulation program (its scheme is presented in Appendix [B]) , containing 
the geometry of WASA detection system and taking into account detection 
and reconstruction efficiencies. The Monte Carlo simulations were carried out 
for the deuteron beam momentum range Pbeam £ (2.127,2.422) GeV/c corre¬ 
sponding to the experimental ramping. The program first chooses randomly the 
momentum magnitude of the nucleon inside the deuteron beam and deuteron 
target, respectively, according to the Fermi distribution m. The direction 
of the momentum vector is chosen isotropically. Then, the total proton-proton 
invariant mass and the vector of the center-of-mass velocity are deter¬ 
mined. Next, the effective proton beam momentum is calculated in the 

frame where one of the proton is at rest and momentum of protons is generated 
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isotropically in the proton-proton center-of-mass frame. Further on, the mo¬ 
menta of outgoing particles are transformed to the laboratory frame and are 
used as an input in the simulation of the detection system response with the 
GEANT computing package. Detailed description of the next simulation steps 
is presented in Appendix]^ For each of A^o simulated event, we assign a weight 
corresponding to the differential cross section, which is uniquely determined 
by the scattering angle and the total proton-proton collision energy 


^pp- 


The factor Nq/2tt in Eq. (8.7) is a normalization constant. It results from 
the fact that the integral is not dimensionless and its units correspond to the 
units of the cross sections used for the calculations. Therefore, it must be 
normalized in such a way that the integral over the full solid angle is equal to 
the total cross section for the elastic scattering averaged over the distribution 
of the total proton-proton invariant mass resulting from the Fermi distri¬ 
bution of the target and beam nucleons. In the absence of the Fermi motion it 
should be simply equal to a total elastic cross section for a given proton beam 
momentum. A factor 27r comes from the fact that protons taking part in the 
scattering are indistinguishable. 

The differential cross section for quasi free dd —>■ ppngpUgp reaction is a 
function of the scattering angle 6* and the total energy in the proton-proton 


centre-of-mass system 


^pp 


which is dependent on effective proton beam mo¬ 


mentum seen from the proton in the proton-proton system. In order 

to calculate it, we have used the cross section values for proton-proton elastic 
scattering pp pp computed based on the SAID program El because the 
EDDA collaboration data base El is insufficient in the region of interest. 
The distribution of the effective beam momentum, as well as a comparison 
of the SAID calculations and the existing differential cross section from the 
EDDA measurements are shown in Fig. 


As we can see, the differential 
cross sections calculated using the SAID programme are in agreement with 
distributions measured by the EDDA collaboration. 


The differential cross section for appropriate and 6* 


was calculated 


using bilinear interpolation in the — 0* plane according to the formula: 


da 

dQ 



(1 - 0(1 - »•■) + f(i - »•■)+ 

+(1 - 


( 8 . 8 ) 


where t and u variables are defined in right panel of Fig. 8.8 
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0 * 



Figure 8.8: (left) Differential cross sections for proton-proton elastic scattering 
as a function of the effective beam momentum f^aln ^ three values of the 
scattering angle 9* in the CM frame. Black points show EDDA collaboration 
while lines denote SAID calculations 


data 120 


119 . Distribution of the 


effective beam momentum for quasi-free pp —?• pp reaction calculated for the 
deuteron beam momentum range p(,earn £ (2.127,2.422) GeV/c is also presented 
in the figure, (right) Bilinear interpolation of the differential cross section 

dO y^be. 


■'beam-’ right figure is adopted from 
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The number of experimental events Aexp was determined based on condi¬ 


tions and cuts described in details in reference [^. The analysis is based on 
the events selected by the hardware trigger for elastic scattering (Sec. 4.5.1) 
and we carried out primary events selection applying condition of exactly one 
charged particle in the Forward Detector (FD) and one particle in the Central 
Detector (CD). 

In Ref. we can find detailed studies of the possible background reaction 
contributions. In case of this analysis the dominating background processes 
are dd dTr+UspUsp, dd d^ptUsp and dd —)■ pPspnUsp, where the subscripts 
sp, b and t denote the spectators, particles from the beam, and from the 
target, respectively. In order to reject events corresponding to the charged 
pions registered in the Central Detector, the cut on the energy deposited in the 
Electromagnetic Calorimeter (SEC) vs. energy deposited in Plastic Scintillator 


Barrel (PSB) spectrum was applied and is presented in Fig. 8.9 


It is not possible to separate quasi-elastic pp scattering from the quasi¬ 
elastic dp scattering, however it was investigated that for the forward scattering 
angles of about 9fd = 17°, the dp cross sections are about 20 times smaller 
than pp cross sections and we take this uncertainty of about 5% as a systematic 
error to the final result. The applied cut in polar angle 9fd is shown in 
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Edep(PSB) [GeV] 


Figure 8.9: Experimental spectrum of the energy loss in the Plastic Scintillator 
Barrel shown as a function of the energy deposited in the Electromagnetic 
Calorimeter. The applied cut is shown as a black line. Points corresponding 
to pions are concentrated for Edep(PSB) around 0.002 GeV. 


Pig. 8.10 In order to subtract the background coming from dd —>■ Pbdtn 


sp 


reaction, the range 6cd G (40,100)° was taken into account in further analysis. 
Additionally, the background was subtracted in Z\0 = (t)FD — 4>CD spectrum. In 
order to symmetrize the background instead of |Z\0| we define (27r+Z\(/))mod27r. 
Afterwards, the background was fitted with 1st order polynomial for each 
of excess energy Q intervals. The exemplary (27r + A0)mod27r spectrum is 


presented in Fig. 8.11 



0fd Id®gl 0FD 


Figure 8.10: Correlations between the polar angles 9fd and 9cd of dd 
ppngpHgp reaction as simulated (left panel) and obtained in experiment (right 
panel). Applied cut is marked with red fine. The indicated area correspond to 
the: a) dd —)■ pprigpUsp, b) dd dbPtUsp and dd —>■ ppngpUgp, c) dd pPspn-Usp, 

d) dd PbdtUsp. 
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Figure 8.11: Distributions of (27r + Z\0)mod27r, where Acf) = (f)cD — 4’fd is 
the difference of azimuthal angles in Central Detector and Forward Detec¬ 
tor. The example spectrum for one of the Q intervals (black line) with marked 
fit function (red line) and signal peak after background subtraction (blue line) 
is presented. 


After all cuts and application of conditions described above, the number 
of experimental events was determined and the luminosity was calculated ac¬ 
cording to formula (8.7) for each excess energy interval. In the calculations 
the prescaling factor of the applied experimental trigger equal to 4000, as well 
as shadowing effect equals 9% were taken into account. The latter results from 
the fact that proton is shadowed by the neutron inside the deuteron which 
reduces the probability of the quasi-elastic scattering. Unfortunately, there 
are no experimental results about the shadowing in dd —>■ ppngpUgp collisions. 
However, we can estimate it based on the probability that a neutron shadows 
the proton in one deuteron which equals 0.045 mi and assume that shadow¬ 
ing appears independently in deuteron beam and deuteron target. The rough 
estimation of the probability that the shadowing will not take place in dd 
reaction (1 - 0.045)^ gives about 0.91. 

The final result is presented in Fig. |8.12[ The statistical uncertainty of each 
point is about 1%. The luminosity variation (increase in the excess energy 
range from about -70 MeV to -40 MeV, and then decrease) is caused by the 
change of the beam-target overlapping during the acceleration cycle and also by 
adiabatic beam size shrinking [114 . The obtained total integrated luminosity 
within its statistical uncertainty is equal to L^dd-^ppnspUsp ~ (1326±2) nb“^. For 
further analysis the luminosity was fitted by third degree polynomial aQ^ + 
bQ^ + cQ + d. The fitted function is marked with the red line in Fig. 8 .12, 
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Figure 8.12: Integrated luminosity calculated for experimental data for quasi- 
free dd ppngpUgp reaction (blue points) with fitted third degree polynomial 
function (red line). 


8.3 Systematics 


In case of the dd —?• ^Hen reaction one source of the systematic error originates 
from the definition of the cuts used for separation of high-energetic helium 
in Forward Detector and is equal to about 2%. Additionally we have also 
taken into acconnt an uncertainty due to the method used for the background 
subtraction amounting to 1.6%. Another source of the luminosity calculation 
error is connected to normalization to SATURNE experiment and originates 
from three independent sources: (i) statistical error of the SATURNE data 
(6.5%), (ii) normalization uncertainty of the SATURNE data for the dd 
^Hen cross sections (7%) and (iii) assumption of linear interpolation between 
SATURNE points used for the estimation of the correction for the parametrized 
cross section presented in Eig. |8.3| (<1.8%). 

The systematic errors for dd —?• ppngpHsp analysis resulting from the cuts 
used for the separation of the quasi-free pp scattering from the backgronnd 


(Eig. |8.9| and Eig. 8.10) is equal to about 4.1%. Another contribution to 
the systematic error comes from the assumption of the potential model of the 
nucleon bound inside the deuteron and is equal to about 0.8%. This uncer¬ 
tainty was established as the difference between results determined using the 


Paris 118 and the CDBonn 122 potentials. The next source of the systematic 
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error may be attached to the assumption of the bilinear approximation of the 
cross section shown in Fig. 8.8 (right). This systematic uncertainty was esti¬ 
mated using assumption in which instead of the interpolation we took the cross 
section value from the closest data point in the effective proton beam momen¬ 
tum - scattering angle plane. The performed calculations give the difference 
of about 1.8%. Additionally we have also taken into account an uncertainty 
related to the background subtraction in (27r -|- A0)mod27r spectra which does 
not exceed 0.6%. The systematic uncertainty includes also contribution con¬ 
nected to the shadowing effect. Until now, we have no theoretical estimation of 
the possible error of this effect, therefore conservatively we take as an system¬ 
atic uncertainty half of this effect: 4.5%. In the systematic error calculation 
we also take into account the uncertainty 5% resulting from the background of 
the quasi-elastic dp scattering. The normalization error includes also normal¬ 
ization uncertainty of the EDDA data (4%) and the systematic errors for pp 
elastic scattering cross-sections (2.7%) mg. The cross section was approxi¬ 
mated by the calculation using the SAID procedure. Because, the SAID cross 
section very well describes EDDA data, we assume the systematic errors of the 
differential cross section based on EDDA calculations. 

The total integrated luminosity determined based on dd —?• ^Hen and the 
quasi-free dd —?• ppngpHgp reactions with statistical, systematic and normal¬ 
ization error are equal to 
and = (1326 ± 2„„ ± 108 

systematic and normalization errors were calculated by adding in quadrature 
the appropriate contributions described above. 

To summarize, the luminosity was calculated based on two reactions: dd 
^Hen and the quasi-free dd —>■ ppUspUsp- The obtained results are consistent, 
however within large normalization errors. 


-‘dd-^^Hen ~ (1102 ± 2stat =■= ^Ogyst 


zb 2Sfiijst i ^^^norm} ^0 


-1 


syst ± Q'^norm) nb \ respectively. The 






Chapter 9 

Results and interpretation 


This chapter presents the determination of the excitation function for dd 
^Hen7r° reaction. Moreover, the contribution from the main experimental back¬ 
ground in the "Signal Rich" and "Signal Poor" region was investigated. The 
obtained results are described in sections below. 


9.1 Excitation function 


The excitation function for dd —> ^Hen7r° process was determined for the region 
where we expect the signal ("Signal Rich" region) after all applied conditions 
and cuts presented in Chapter The excitation curve is shown in left panel 
of Fig. |9.1[ 




Figure 9.1: (left) Excitation function for the dd —> ^Henyr® reaction including 
events selected via applied conditions and cuts, (right) Excitation function ob¬ 
tained by normalizing the events selected in individual excess energy intervals 
by the corresponding integrated luminosities and efficiencies. 


The number of events in each excess energy interval was divided by the corre¬ 
sponding integrated luminosity L{Q) determined based on quasi-free proton- 
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Results and interpretation 


proton scattering (Sec. 8.2) and corrected for total efficiency (Chapter [^. The 
normalized excitation function is presented in the right panel of Fig. |9.1[ 


9.2 Upper limit of the total cross section 


The shape of obtained excitation function can be well described with a quadratic 
function fit resulting in the value per degree of freedom of 1.1 (Fig. 9.2). 
The excitation function does not indicate significant sharp enhancement for 
energies below the r] production threshold which could be interpreted as a 
resonance-like structure. Therefore, we can only determine an upper limit for 
the cross-section for formation of the ^He-77 bound state and its decay into the 
^Hen7r° channel. The excitation function was fitted with quadratic function 
describing the background combined with the Breit-Wigner function which can 
account for the signal from the bound state: 


a{Q,r,Bs,A) 

where: Bg -binding energy, 
r - width, 

A - amplitude. 


A ■ rV4 

(Q-Bgf + r^/4’ 


(9.1) 


In applied fit the polynomial coefficients and the amplitude A of the Breit- 
Wigner distribution are treated as free parameters while the binding energy 
Bg and width F are fixed parameters. The fit was performed for various values 
of Bg and F. The binding energy and the width were varied in the range from 
0 to 40 MeV and from 5 to 50 MeV, respectively. An example of the fit for 
T = 40 MeV and Bg = 30 MeV is presented in Fig. |9.2[ 

The upper limit of the total cross section was calculated as: 


upp 

'CL=90% 


= k-aA, 


(9.2) 


where a a is uncertainty of the amplitude A obtained from the fit and k is 
statistical factor equal to 1.64485 corresponding to the confidence level (CL) 
of 90%. The example values of the obtained upper limit are given in Table [9Tj 
As one can see, the upper limit depends mainly on the width of the bound 
state while its dependence on the binding energy is only slight. The obtained 
upper limit as a function of bound state width is presented in Fig. |9.3| for 
binding energy 30 MeV. Its value varies between 21 to 36 nb. The green area 
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denotes the systematic errors which contributions are described in details in 
the next section. 



Figure 9.2: Excitation function for the dd —>■ ^Hen7r° reaction obtained by nor¬ 
malizing the events selected in individual excess energy intervals by the cor¬ 
responding integrated luminosities and corrected for acceptance and effciency. 
The red solid line represents a fit with second order polynomial combined with 
a Breit-Wigner function with fixed binding energy and width equal to 30 and 
40 MeV, respectively. The blue dotted line shows the second order polynomial 
corresponding to the background. 


Bg [MeV] 

r [MeV] 

[-b] 

Bg [MeV] 

r [MeV] 

[-b] 

10 

10 

22.05 

30 

10 

21.04 

10 

20 

21.34 

30 

20 

21.35 

10 

30 

24.32 

30 

30 

24.32 

10 

40 

29.26 

30 

40 

29.27 

10 

50 

36.13 

30 

50 

36.15 

20 

10 

22.61 

40 

10 

21.55 

20 

20 

22.72 

40 

20 

20.08 

20 

30 

27.00 

40 

30 

22.12 

20 

40 

33.93 

40 

40 

25.81 

20 

50 

39.24 

40 

50 

30.94 


Table 9.1: The upper limit of the total cross-section for the dd (^He- 
v)bound ^Henvr^ process determined at CL=90% for different values of bind¬ 
ing energy Bg and width B. 
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Results and interpretation 



r [MeV] 


Figure 9.3: Upper limit of the total cross-section for dd (^He- 77 )bound ^ 
^Hen7r° reaction as a function of the width of the bound state. The binding 
energy was set to 30 MeV. The green area denotes the systematic uncertainties. 


9.3 Systematics 


Systematic studies were carried out based on Ref. [^. The variation of the 
selection conditions by ±10% gives the systematic error which equals about 
6% (the highest contribution to the error comes from changing the range of the 


"Signal Rich" region A (see Sec. 6.1.3). Another, significant contribution to 
the systematic error of the upper limit is connected to the luminosity determi¬ 


nation based on quasi-free pp reaction (Chapter and equals about 8% and 
5% which correspond to the systematic and normalization error, respectively. 
Additional source of systematic errors comes from the assumption of the Fermi 
momentum distribution of nucleons inside the ^He. The distribution was used 
in Monte Carlo simulation of the bound state production and decay (Chap¬ 
ter [^. Current analysis was performed with the Fermi momentum distribution 
based on AV18-TM potential model in Another available momentum dis¬ 
tribution based on CDB2000-UIX model only slightly changes the acceptance 


for simultaneous registration of all particles in WASA detector (see Sec. 5.3) 
providing systematic error of about 1%. The uncertainty caused by the fit of 
quadratic or linear function to the background, estimated as changes 

from about 5% (T=5 MeV) to 25% (±=50 MeV). The systematic error is cal¬ 
culated by adding in quadrature all contributions described above and varies 
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from 12% to 27% as shown by green area in Fig. 


9.3 


9.4 Background studies 

The excitation function presented in the previous section is dominated by 
the background which interferes in the dd —?• hound ^ ^Hen7r° process. 

Therefore, understanding of the background processes is crucial for the consid¬ 
ered investigations. We performed the studies of two reactions being the main 
background contributions: 

1. dd —> ^Hen7r°, 

2. dd ^ ^HenyrO. 


The simulations of above reactions were carried out for the beam mo¬ 
mentum range corresponding to the experimental range Pbeam G(2.127,2.422) 
GeV/c. The first reaction was simulated according to direct production with 
the uniform distribution over the phase space. The second process proceeds 
with excitation of an N* resonance which subsequently decays in the nucleon- 
7r° pair. The detailed description of the simulations of dd —>■ ^HeA^* —> ^Hen7r° 
reaction can be found in the Appendix 

We compared the excitation curves obtained from simulations of back¬ 
ground processes with the excitation function determined from the experi¬ 
ment. The excitation function for experimental data normalized over luminos¬ 
ity was fitted with a A ■ WMCbcgi + B ■ WMCbcg 2 function where WMCbcgi 
and WMCbcg 2 are the excitation functions obtained from the simulations of 
dd —7- ^Henvr® and dd ^HeA”* ^Henvr® reactions, respectively. The simu¬ 
lations were carried out using the conditions and cuts described in Chapter 
The result of the fit is presented in Fig. |9.4| for region where we expect the sig¬ 
nal from the bound state G (0.07,0.2) GeV/c) and for the region poor in 
signal (pajjg G (0.3,0.4) GeV/c). It is important to stress that the experimental 
excitation function for the "Signal Rich" region cannot be well described only 
by the background contributions (see left panel of Fig. 9.4). It can indicate 
missing contribution from another background processes, or the influence of 
wide ^He-r 7 or ^He-A* bound state. In contrast, as shown in the right panel of 
Fig. |9.4[ the experimental excitation function from the "Signal Poor" region 
is very well described by the background reactions. 
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Figure 9.4: Experimental excitation functions (red circles) fitted with two 
background reactions: dd ^Hen7r° (green squares) and dd —>■ ^HeA^* 
^Hen7r° (magenta squares). A sum of both background contributions is shown 
as blue triangles. Left and right panels show results for the "Signal Rich" and 
"Signal Poor" regions, respectively. 













Chapter 10 

Conclusions and outlook 


This dissertation describes the search for the ^He-ry bound state via the study 
of the excitation function for the dd —> ^Hen7r° reaction. It includes detailed 
description of the performed experiment, the data analysis and obtained re¬ 
sults. 

The measurement of dd ^Henvr® reaction was performed in 2010 with 
the WASA-at-COSY detection system using the deuteron beam and deuteron 
pellet target. To reduce the systematic uncertainties of the beam momentum 
value, the ramped beam technique was used. The experiment was performed 
for the beam momentum range from 2.127 GeV/c to 2.422 GeV/c correspond¬ 
ing to the excess energy range of Q G (-70,30) MeV for the ^He-ry system. 

The performed analysis allowed for the determination of the excitation 
function for dd —)• ^Hen7r° process and the estimation of the upper limit of the 
cross section for the ?y-mesic ^He formation and decay. Events corresponding 
to production of the ry-mesic bound state were selected via cuts on the ^He 
momentum, neutron and 7r° kinetic energies as well as on the opening angle 
between n-7r° in the center of mass frame. The cuts were adjusted based on 
Monte Carlo simulations of the {"^He-r])bound production and decay. The simu¬ 
lations assume that the ry meson inside bound state is absorbed by the one of 
neutron and excites it to N* resonance which subsequently decays into n-7r° 
pair. The ^He plays here the role of spectator moving with Fermi momentum 
in the center of mass system. 

The total integrated luminosity in the experiment was determined based 
on the dd —)• ^Hen and quasi-free pp —)• pp reactions. It amonnts to ~ 

(1102±2stof±285y5f±107„o,rm) nb“^ and = (1326±2stot±1085y5f± 

Q4:norm) nb“^, respectively. The obtained results are consistent within system- 
atics and normalization error. The excess energy dependence of luminosity 
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Conclusions and outlook 


used for relative normalization of the dd —?• ^Henvr® reaction events was deter¬ 
mined based on quasi-free pp pp reaction for which the WASA acceptance 
is a smooth function of the beam momentum. 

The obtained excitation function of the dd —>■ ^Hen 7 r° reaction does not re¬ 
veal the resonance-like structure which could be interpreted as the indication 
of the ^He-r] bound state with width less than 50 MeV. The fit of the exci¬ 
tation curve with the Breit-Wigner distribution combined with second order 
polynomial allows to determine the upper limit of the dd bound 

^Hen 7 r° cross section as function of the bound state width and binding energy. 
The upper limit varies from 21 to 36 nb for the width varying from 5 MeV to 
50 MeV. The obtained upper limit is by factor of five larger than the theoreti¬ 
cally estimated value of the cross section for the dd —>• {"^He-r])bound ^Hepyr" 

reaction |^ . Therefore, we can conclude, that the current measurement does 
not exclude the existence of ^He-r 7 bound state in considered process. 

The search for r]-mesic ^He is complex due to huge background contribu¬ 
tion. This background largely comprises two processes dd —?• ^Henvr® and 
dd —>■ ^HeN* ^Hen 7 r° being direct production and production via N* res¬ 
onance without formation of {'^He-r])bound, respectively. The contribution of 
considered processes to the experimental excitation function was investigated 
for the regions rich in signal ("Signal Rich") and poor in signal ("Signal Poor") 
selected in the ^He momentum. In case of the region rich in signal the com¬ 
bination of the considered background reactions is not sufficient to properly 
describe the experimental data. This may suggest an influence of some other 
background process not taken into account, or the sign of the very wide ^He-r] 
or ^He-V* bound state. This result is a subject of interpretation of few theo¬ 
retical groups miiiii- 

In 2014 the search for 77 -mesic bound states at WASA was extended to ^He-r^ 
system. Based on the new research hypothesis about the mechanism of the 
decay of r^-mesic nucleus we have elaborated an experimental proposal m 
for the search of the ^He -77 bound state which was accepted by the COSY 
Advisory Committee in February 2014 and already in the year 2014 we com¬ 
pleted successfully the experimental run. Analysis of the collected data is in 
progress. The search for p and p' -mesic nuclei is carried out also by another 
experimental groups, e.g. at J-PARC and at GSI |125 , 126j . 













Appendix A 
The 77 meson 


The rj meson was discovered in 60’ by Pevsner m as a resonance in three 
pion invariant mass spectrum. Its properties were investigated for many years 
and are summarized by the Particle Data Group |128| . The main information 
about the 77 meson is briefly presented in Table |A. 1 [ 


mass 

547.853 ± 0.024 MeV 

width 

1.30 ± 0.07 keV 

jG(jPC) 

o+(o-+) 

Decay modes 

Branching ratio 

Charged modes 

28.10±0.34 % 

T] Tr~^Tr~7r^ 

22.74±0.28 % 

7] TT~^TT~^ 

4.60±0.16 % 

other modes 

0.76 % 

Neutral modes 

71.91±0.34 % 

77 ^ 27 

39.31±0.20 % 

77 ^ 37 r° 

32.57±0.23 % 

other modes 

0.03 % 


Table A.l: rj meson properties and its main decay modes. The data are taken 
from |128| . 
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The 7 ] meson 


77 is a neutral meson with zero spin J and negative parity P. Together with 
77 ', pions: 7 r°, tt"'', 7 r“ and kaons: K^, , K~^, K~ it is a part of pseudoscalar 

meson nonet which is schematically presented in Fig. The mesons are 
arranged here according to strangeness S along the Y axis and according to 
the isospin component I 3 along the X axis. 



Figure A.l: Multiplet of pseudoscalar mesons. The horizontal axis corresponds 
to 3’’'^ component of the isospin p while the vertical axis corresponds to 
strangeness S. 

In the Gell-Mann quark model mesons consist of quark-antiquark pairs. 
According to SU(3) flavour symmetry of three lightest quarks u, d and s the 
mesons form a singlet and an octet with the following quark contributions: 

771 = ;^ {dd + uu + ss ), 

778 = -^{dd -\- uu — 2 ss). 

The observed 77 particle is the combination of the 771 and 773 states: 


I 77 >= rjgcosO — rjisinO, 


where 6 is mixing angle determined experimentally and equals 0 =-15.5°±1.3° 
According to Bass and Thomas the flavour singlet component 771 can mix with 


pure gluonic states 16,17 what is important from the point of view of the 77 
















85 


meson properties inside the nuclear matter. 

Many strong or electromagnetic decay channels of rj meson is forbidden 
by C, P, CP or G symmetry conservations. The decay into two or four pi¬ 
ous is forbidden due to P and CP invariance while charge conjugation does 
not allow to occur the r] 7 r° 7 , r] 7 r° 7 r °7 and rj —?• 7 r° 7 r° 7 r °7 processes. 
G-parity conservation forbids the decays into three pions, however they occur 
with isospin violation and are dominant processes together with second-order 
electromagnetic —>■ 77 decay (Table A.l). The investigation of rare rj sym¬ 
metry violating decay processes is very important for precise studies of the 
QCD symmetries. 




Appendix B 

Simulation of dd —^ ppnsptigp 
reaction 


The quasi-free pp reaction described in Sec. |8.2| is schematically presented in 
Fig. |B.l|and its simulation is described below. 


before the reaction 


after the reaction 



deuteron beam deuteron target 


quasi-free protons 

m 


neutron spectators 


Figure B.l: Schematic picture of the quasi-free pp —?• pp reaction. Red and blue 
circles represent protons and neutrons respectively. The Fermi momentum of 
the nucleons inside the deuteron beam and deuteron target is presented by the 
dotted arrows and the deuteron beam momentum by the dashed one. 


1- Pbeani IS generated with uniform probability density distribution in the 
range oi pbeam G (2.127,2.422) GeV/c. 
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2. The deuteron beam (6) as well as the deuteron target (t) are considered 
as a proton-neutron bound systems {pb,nb) and (pt,nt), respectively. The 
neutron nb and Ut momentum vectors are distributed isotropically in 
the spherical coordinates of the deuteron beam and the deuteron target 
0Ff,t) with Fermi momentum distribution of nucleons inside 
deuteron Hzl and transformed into Cartesian coordinates. Fermi mo¬ 
mentum distributions of proton and neutron bound inside a deuteron 
derived from two different potential models, namely PARIS |118| and 
CDBonn |122| are shown in Fig. B.2, The neutrons four-momenta are 
calculated based on spectator model assumption (llPn^tP = ^^e 

center of mass (CM) frame and transformed using Lorentz transforma¬ 
tion into the laboratory frame. 



Figure B.2: Fermi momentum distribution of nucleons inside the deuteron for 
PARIS (full line |118| and CDBonn (dashed line) |122] potentials. The distri¬ 
butions were normalized to unity in the momentum range from 0 to 300 MeV/c. 


3. Based on Fermi momentum pp^ ^ values, proton Pb mass in the deuteron 
beam and proton pt mass in the deuteron target are calculated according 
to equation (|B.1|). 


^Pb,t = + ^l- 2md^m2 + ' . (B.l) 

Four-momenta of both of protons are calculated in the CM systems of 
deuterons and are transformed into the laboratory (LAB) frame. 



















Simulation of dd —> ppUspUsp reaction 


4. The proton-proton invariant mass is calculated in LAB system ac¬ 



cording to equation (B.2): 



(B.2) 


where Ep ^, Ep^ and Pp ^, Pp^ are energies and momenta in LAB frame for 
Ph and pt, respectively. 

5. Based on Pb momentum is calculated in the system where pt is at 

rest: 



(B.3) 


6 . Four-vectors of protons pb and pt are transformed into the proton-proton 
CM frame. 

7. Scattering between protons is considered. Scattering angle 6 * as well as 
azimutal angle (p* are simulated isotropically. Four-momenta of scattered 
protons are calculated and transformed to the laboratory system. 


8 . WASA Monte Carlo (simulation of the detection system response by 
GEANT package) is carried out for generated events. 






Appendix C 
Simulation of 

dd —> ^HeiV* —> ^HenTT® reaction 


The simulation of dd —>■ ^HeA* —>■ ^Henyr^ reaction, being one of the main 
processes contributing to the background (Sec. |9.4 ), was carried out for the 
beam momentum range Pbeam £ (2.127,2.422) GeV/c corresponding to the 
experimental ramping. The main assumptions of simulation are schematically 
described in following points: 


1. The deuteron beam momentum value Pbeam is generated with uniform 
probability density distribution in range of Pbeam G (2.127,2.422) GeV/c 
and then the square of invariant mass of the whole system Sdd is calcu¬ 


lated using Eq. (5.2) presented in Sec. 5.2 


2 . The invariant mass ^/sM is distributed randomly according to the distri¬ 
bution presented as follows: 


rvVmax 

= / PS{W) ■ BW (V^ - IT - msHe, At*, • dW, 

^ yVmin 

(C.l) 

where: 

•W = - mjv* - msHe 


is the excess energy avaliable in the CM frame with minimum and max¬ 
imum values equal to kTnm = 0 and Wmax = -s/^dd — — 'ITt-n — 

respectively; 
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Simulation of dd —)■ > ^HeriTr” reaction 


• PS{W) = + msHe] [sm - {m^* - m3^e)f 

= [2Vsdd - 1^] [Sdd - (v^ -W- 2msHe ?]/ (27s^3) 


is a Phase Space factor for 2-body reactions which is proportional to '/W 
near the rj production threshold and to 1/W above the threshold CT; 

. BW (ySS -W- mme, r„.,E^.) = = 

__ 

-Sjv* /4 


is a Breit Wigner distribution of N* resonance with energy £'jv*=1535 MeV 
and width P/v* = 150 MeV. The BW distribution is presented schemati¬ 


cally in Fig. C.l while cr(y^s^) distribution in Fig. C.2 



Figure C.l: Breit Wigner distribution of N* resonance with energy 
F'Ar*=1535 MeV and width /V*=150 MeV. The green arrows show the 
maximum resonance masses mjv* and mj^* for the beam momentum 
Pbeam = 2.127 GeV/c and Pfteam = 2.422 GeV/c, respectively. Blue arrow shows 
the sum of pion and neutron masses {m -,^0 and m„) which is the lower limit of 
the resonance mass in dd —>■ ^HeV* —>■ ^Hen7r° process. 

3. Excess energy available in the CM frame W is distributed according to 
the PS{W) ■ BW — W — msHei -Cv*, Ej^*) distribution. 
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Figure C.2: a{^/sM) distribution determined for the considered beam momen¬ 
tum range. The dashed line denotes the minimal and maximal values of Pbeam 
which are equal 2.127 GeV/c and 2.422 GeV/c, respectively. 


4. The resonance mass m^* is calculated as m^* = ^/sm — W — msHe and 
is limited, because of two conditions: 

• mjsi* + msjje < (the whole available energy is used to produce N* 
and ^He), 

• mjsi* > m^o -h rrin (resonance mass should be enough to decay into 
neutron and 7r°). 


5. The neutron and pion momentum vectors are simulated isotropically in 
the N* frame in spherical coordinates and transformed into Gartesian 
coordinates. The absolute value of neutron and pion momenta is 


hxed by equation (5.5) described in Sec. 5.2 


6 . The gamma quanta are simulated isotropically in the 7r° frame in spher¬ 
ical coordinates with momenta p*** = mT^o/2. 


7. The four-momentum vectors of ^He, neutron and gamma quanta are 
transformed into the center of mass frame and next into laboratory frame 
by means of Lorentz transformation. 


8 . The WMG simulations of the WASA detector response are carried out. 








List of Abbreviations 


WASA Wide Angle Shower Apparatus 

COSY Cooler Synchrotron 

PDG Particle Data Group 

FD Forward Detector 

CD Central Detector 

FWC Forward Window Counter 

FPC Forward Proportional Chamber 

FTH Forward Trigger Hodoscope 

FRH Forward Range Hodoscope 

FRI Forward Range Intermediate Hodoscope 

FRA Forward Absorber 

FVH Forward Veto Hodoscope 

MDC Mini Drift Chamber 

PSB Plastic Scintillator Barrel 

SCS Superconducting Solenoid 

SEC Scintillator Electromagnetic Calorimeter 

WMC Wasa Monte Carlo 

CM Center of Mass System 

LAB Laboratory System 

DAQ DATA Acquisition 

FPGA Field Programmable Gate Array 

QCD Quantum Chromodynamics 
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